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Abstract. This paper models the classical diffusion of a main particle through 
a heatbath by means of a pre-limit microscopic representation of its drifted 
momentum and energy transfers at collision times. The collision point lin- 
ear interpolated path can be approximated by the solution to the "inscribed" 
continuous stochastic differential equation using the same drift function. Em- 
ploying results from stochastic mechanics it is then shown that the combined 
main particle/heatbath system does not exchange or radiate energy if the 
probability distribution for the position of the main particle is derived from 
Schrodinger's equation. Furthermore it is shown that the main particle dis- 
tance traveled between collisions and the mean inter-collision time must satisfy 
a type of Minkowski invariant. Hence if there is a correlation between the pre- 
and post-collision velocities of the main particle through a collision point then 
the mean distance traveled can be related to the mean inter-particle collision 
times via a Lorentz transformation. The last Section shows that this approach 
can be applied to all elastic main particle/heatbath particle collisions either 
via direct calculation involving modeling the collision scattering or by altering 
the properties of the heatbath. 



Introduction 

Ever since Einstein's introduction of the molecular-kinetic theory of heat in 1905 
the Brownian motion/Mar kovian formalism has been applied to a large variety of 
topics including classical particle diffusion and stochastic mechanics. The first sub- 
ject is more the domain of statistical mechanics and focusses on thermodynamic 
properties, Goldstein Isothermal Flows Garbaczewksi j2j, transport equations 
(e.g. Master equation, Boltzmann's equation or Kramer's equation) and Markov 
Chains, for instance Posilicano van Kampen or Gamba The second 
topic falls under the interpretation of quantum mechanics see for instance the re- 
cent review by Carlen [S] or Nelson [?]• For a historical view on the development 
of Brownian Motion consult Nelson The present paper represents a (larger) 
particle (the " main" particle) diffusing through a heatbath of smaller particles em- 
ploying a "finite-energy" Markovian difference equation and investigates various 
energy conditions using results that have been developed in stochastic mechanics. 
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The simplifying assumption for diffusion is that the frequent energy exchanges 
between the main and heatbath particles induce a continual acceleration and de- 
celeration that make the motion of the main particle look macroscopically as if it 
has no memory of its previous whereabouts. Mathematically the Brownian motion 
process is the result of a limiting procedure involving an infinite amount of collisions 
exchanging an infinite amount of energy. The original derivation of the associated 
diffusion equation for the density distribution is due to Einstein jHI- As a result 
of the limiting process the resulting Markovian process is not differentiable and 
therefore no energy or momentum can be ascribed to the main particle. Brownian 
motion is very attractive because of its " random step" intuition and because of the 
tractability of the associated diffusion equation. 

Many ways of associating momentum and energy with the motion of a diffusing 
particle have been suggested in the literature. A straightforward and transparent 
approach is to augment the position process for the diffusing particle with a momen- 
tum process and investigate the position process in the limit that the momentum 
becomes extremely large. Modeling both position and momentum Nelson 8^ showed 
that the position of the main particle converges in probability to the solution of 
a stochastic differential equation if the associated momentum process contains an 
ever larger mean reversion (and variance). The more typical - and mathematical - 
construction of Brownian motion is based on successive probabilistic additions of 
orthogonal functions, see Rogers & Williams iTJ or Karatzas & Shreve jjl]. These 
two constructions do not consider the underlying particle collision process or the 
implied energy exchange between the main particle and heatbath environment. 

This paper uses a direct approach by creating a pre- limit microscopic representa- 
tion of a main particle by modeling its momentum and energy transfers at collision 
times. The assumption is that the collisions occur at stopping times tj, j > with 
tj+i — tj = Tj, j > 0, so that the main particle will have corresponding positions 
Xj = X {tj,l3) G R", j > 0, /? = 2/r at these collision times. Here (3 denotes the 
main particle path dependence on the mean particle inter-collision time t. The 
particle position at the collision point Xj+i depends on the position of the previous 
collision point Xj, a drift term {xj, tj) and a Gaussian random shock. Since the 
(main) particle drift over the interval [tj, ij+i] , j > 0, is a function of Xj, tj,j > 0, 
the position process of the main particle Xj,j > 0, reduces to a discrete Markovian 
stochastic process. The crucial assumption in this paper is that the inter-collision 
times TjJ > 0, can be represented by a second (or higher) order gamma distribu- 
tion. The main reason for this assumption is that the energy and momentum of the 
main particle between collisions are then well defined and finite with probability 
one. 

The set of coUision points {xj = x {tj,f}) |j = 0, 1, ...} does not prescribe the po- 
sition of the main particle at an arbitrary time t ^ {tj\j = 0, 1, ...}. A reasonable 
estimate for the position of particle at time t is to define x {t, /3) ,t > 0, as the lin- 
ear interpolation on the collision positions. So x (t, (3), is linearly interpolated from 
x{tj,l3) and x(tj+i,/?) where tj <t < tj+i (to — 0, the origin). The expectation 
is that if the number of collisions increases the process x{t,(3) S K", approaches 
the continuous stochastic process x (t) G K", t > 0, which satisfies an appropriate 
stochastic differential equation with drift b'^^x^t). This obviously depends on the 
properties of the drift function, the variance function in the stochastic shock and 
the existence of collision points. Section 1 introduces a set of conditions for which 
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the discrete process x{t,f3) converges a.e. to the strong solution of a stochastic 
differential equation. 

By design, if the particle experiences a collision at time t with x{t, (3) = x then 
the next collision will occur at time t+T2 while the last collision of the main particle 
occurred at time t — ti. Hence the main particle has a (forward) momentum (and 
energy) towards the next collision equal V2 = {x{t + t^) — x) /t2- The forward 
momentum is a random variable depending only on the current position of the 
main particle due to the fact that the stochastic process x{tj),j > 0, is Markovian. 
Similarly, the main particle will have a backward (or incoming) momentum (and 
energy) from the last collision equal vi — {x — x{t ~ ti)) /ti since t — ti is the 
collision time previous to t. The backward momentum is a random variable that 
will have to be conditioned on the fact that it is now known that x{t) — x. At least 
in principle, the distribution for x(t — ti) follows from the distribution of the main 
particle for x{t) employing Bayes' theorem. The forward and backward time step 
perspective for continuous stochastic processes can be gleaned from the extensive 
work on stochastic processes by Nelson Carlen ^21: QHIi Guerra ^1], ^HI and 
see the references in Carlen 0. 

Now that pre- and post-collision momenta and energies are determined for the 
main particle the collision process and energy exchange can be investigated for the 
main and heatbath particle. Section 2 investigates the consequences of an elastic 
collision and introduces the canonical linear relationship between the pre- and post 
collision velocities of the main and heatbath particle. This linear relationship is pa- 
rameterized by a random anti-symmetric matrix Z (manufactured from a random 
unitary matrix U) which incorporates the random center of mass line and collision 
impact angle information. This matrix will be referred to as the collision scatter- 
ing matrix and collision energy exchanges with Z = are referred to as 'simple' 
collisions. Notice that in one dimension all collisions are simple. 

Section 2 shows that the combined kinetic energy Tik of the main and heatbath 
particle for simple elastic collisions can be expressed in the form of a simple qua- 
dratic combination of the forward and backward momenta/velocities of the main 
particle. In one dimension this relationship is exact while in more than one dimen- 
sion the total kinetic energy Hk of the colliding system refers to the energy of the 
motion along the center of mass line of the collision. The remaining kinetic energy 
of the main and heatbath particles is embedded in motion perpendicular to the 
collision center of mass line and remains invariant under the collision. Section 2 
only investigates the canonical case where TCk = "Hk {Z = 0) while Section 4 shows 
that the case TCk = Ti-k (Z) , Z ^ 0, can be reduced to the simple collision case. 

Using this result Section 2 then shows that if the total kinetic energy of the 
system is not conserved then either the main particle is radiating energy into the 
heatbath or the main particle is absorbing energy from its heatbath surroundings. 
Such an exchange is possible as the result of external forces in the form of a potential 
<i>p but in that case the combined kinetic energy Tik and $p together must be a 
conserved quantity, i.e. its expected value over all paths and positions must be 
a constant in time. The main result of the paper is that the only probability 
density that renders the total energy E [Hk + $p] time invariant is the squared 
modulus of the wave function obtained from Schrodingers equation. Appropriate 
conditions and examples for the potential will presented in Sections 2 and 4. This is 
a purely classical representation and Planck's constant is now replaced by a constant 
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depending on the variance of the underlying stochastic process and the main 
particlc/hcatbath particle mass ratio 7. This Section also shows (using collision 
elasticity) that the backward and forward momenta for the heatbath particle are 
tightly correlated if the main particle follows a Markovian path. 

The main particle/heatbath particle elastic collision representation also provides 
a similar quadratic expression showing that the forward and backward velocities of 
the main particle are directly related to the (forward and backward) velocities of 
the heatbath particle. This relationship will often be referred to as the "momen- 
tum" constraint. Section 3 will show that if the heatbath particles are in energetic 
equilibrium with the main particle then the drift of the main particle and the 
correlation between the forward and backward velocities of the colliding heatbath 
particles must depend on the average inter- collision time. The relationship between 
the mean inter-particle collision time and the particle energy can be expressed in 
the form of the geometric Minkovski invariant. A consequence is that the mean 
distance traveled for the main particle and its mean inter-particle collision time can 
then be related via a Lorentz transformation. 

The last Section focusscs on the non-simple elastic collisions and qualifies the 
main and heatbath particle total energy dependence on the anti-symmetric scatter- 
ing matrix Z = Z{U). The same conservation of total energy employed in Section 2 
now produces an equation which also contains the dynamics of the Z matrix. The 
conservation result is discussed and some examples presented but no full solution 
for this case can be derived. However, the final results in this Section show that 
the total energy conservation and "momentum" constraint can be reformulated in 
terms of a transformed heatbath with heatbath particles that have a higher kinetic 
energy and a different correlation structure. It is possible therefore to transform 
the scattering matrix Z away by subsuming it into the heatbath. Conveniently then 
all the results of Section 2 and 3 become valid again for the transformed heatbath. 

The work is organized as follows. Section 1 introduces the details of the main 
particle collision representation and shows that the collision point linear interpo- 
late converges a.e. to the solution of a continuous stochastic differential equation. 
Section 2 employs the fact that the particle collisions are elastic to show that a 
non-radiation condition demands that the only acceptable probability density for 
the position of the main particle is derived from Schrodinger's equation. Section 3 
uses the "momentum" constraint to show that the mean inter-particle collision time 
and inter-particle distance traveled satisfy a geometric Minkowski invariant. The 
last Section shows that the results of Sections 2 and 3 can be extended to almost 
all types of elastic collisions. All proofs have been delegated to the Appendices to 
make the paper more readable. 

The notation employed in this this paper uses E[.] or E[.\.\ for expectation 
or conditional expectation respectively. Typically E[.] indicates an expectation 
over all variables between the brackets including the random scattering matrix Z 
present and the forward and backward velocities, also b'^{x,t) : M" x [0,oo) 1 — > 
M" are referred to as the forward and backward instantaneous drifts of the main 
particle and (7(.x, t) : K" x [0, oc') 1 — > M" x E" are the corresponding variance terms. 
Moreover r refers to the inter-particle collision time, t = E[t] and /3 = 2/r. Often 
the reference to the time t implicitly assumes that this time point is a collision time 
for the main particle. Clearly for finite mean collision times a distinction must 
be made between the collision points t G {tj,j = 0,1, ...} and the remaining time 
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points but the paper docs not always complete the analysis. The norm 
the usual distance norm in M". 



indicates 



1. Diffusion and Energy 

This Section introduces the details of the Xj — x{tj,(3), j — 0,1,..., collision 
representation, the converged process x(t), t > 0, and the distribution for the inter- 
particle collision times Tj,j > 0. It is shown that limp^oo x{t, (3) exists under 
certain smoothness conditions and this Section also investigates the correlation 
between the forward and backward main particle velocities. Some results from 
Nelson and Carlcn will be quoted without proof. 

Following the Introduction let the collision points for the main particle be x(tj) e 
M", j — 0, 1, at collision time t e {tj,j = 0, 1, ...}, and assume that the next 
collision will occur in x{t + T2) at stopping time t + T2. The collision previous 
to time t occurred at t — ti when the main particle was in position x{t — ti). 
The collisions must be assumed "real" in the sense that the main particle actually 
interchanges energy with a colliding heatbath particle. The inter-collision times are 
modeled as independent random variables distributed with a second order Gamma 
distribution so that /^(t) = /3^te~'^*, t ^ 0. As a result the time to the next collision 
T2 and the time since the previous collision ti has the following moments 



E[r2] = E[r,] = f 
Jo 

var{T2) = var{Ti) - 



Ph^e-^'dt = - - r. 
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Most importantly for this distribution it is also true that 
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Higher order Gamma distributions could have been employed for the inter-collision 
times or in fact any distribution such that E < 00 would have been suitable 

but there seems little fundamental difference in the analysis. Interestingly, the 
ubiquitous exponential distribution is excluded due to the last restriction but notice 
that the second order Gamma distribution is the distribution of the sum of two 
exponential random variables. 

By assumption above if t is a collision time with main particle position x{t, (3) = 
x G M" then the previous collision occurred at i — ti . In this case ri is a random 
variable conditional on the fact that a collision occurred in the future t with the 
particle in position x{t, f3) = x G K". Unfortunately, the distribution for ri condi- 
tional on this event is no longer a second order Gamma distribution, see Feller [T^ 
for a more complete discussion. The future conditioning alters the distribution of 
Tl which is easy to see since any collision time t — ti previous to t must satisfy 
^ t — Tl ^ t. Hence ^ ri ^ t counter to the domain definition of a Gamma 
distribution. Moreover, there is always the possibihty that no previous collision 
occurred so the real density for ri is in fact defective. However, ioi t :§> t it 
is reasonable to assume that ti is at least approximately a second order Gamma 
distribution with E [n] = E [T2] and E [rf ^] = E [ra"^] . 
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The forward and backward step for the position process x{t, (3) at colhsion time 
t are now defined as foUows 

x{t + T2,I3)- x{t, (3) = A+x(t, 13) = b+{x{t),l3, t)r2 + ctA+z, (1.2a) 
x{t,l3) - x{t~Ti,/3) = A-x{t,l3) = b~{x{t,f3),t)Ti +ctA"z, (1.2b) 

where ti, T2 are the inter-coUision stopping times with (independent) second order 
Gamma distributions with 6± (a;, i) : x [0, oo) i — > R" and cr(x, i) : R" x [0, oo) i — > 
M" X M" as the drift vector and variance matrix respectively. Here P = 2/t — 
^ ['''2"^] ~ ^ ['''r^]' equation p. II) . and A+z = z(tj-|_i) — 2(t) where z{t) £ 
M",j ^ 0, is a Gaussian process and A^z is a Gaussian increment independent of 



A+z with 



Due to the aforementioned issue with the distribution for ti the postulated form 
for the backward velocity in lll.2b|) cannot be correct. For fixed T2 and ti Bayes' 
Theorem implies that A^x{t,(3) must have a Gaussian distribution but not nec- 
essarily with the same variance matrix a = a{x,t,P). As ti and T2 are random 
variables the forward increment is no longer Gaussian and the distribution of the 
backward velocity distribution is unclear. Therefore equation ljl.2b(l should have 
been written as 

x{t,P)-x{t-Ti,l3) ^ A-x{t,(3) ^b-{x{t,f3),t)Ti + o-^A^z, (1.2c) 

with 

E[A-x{t,m<t,P)] 
b^{x(t,f3),t)- , ^^^^^^ 

apA-z = A-x{t,p) - b^{x{t,l3),t)Ti, 

where again i?[A~2] ~ ti. Furthermore for finite /3 the increment A~ z is not 
a Gaussian increment. However the expectation is that for progressively smaller 
average inter-particle collision times (r — > 0) that b'^{x,t) b~{x,t), cr^ — > ct 
and A~z becomes approximately Gaussian. In other words bp{x,t) approaches the 
backward drift used in stochastic mechanics and for sufficiently large /3 (small r) 
equation (jl.2dp approaches ljl.2b|l . Equation H1.2a|l is therefore a definition while 
H1.2b|) is only true in the limit of sufficiently small time steps. 

Due to H1.2a|l . (jl.2b|) the forward and backward velocities can now be written as 

V2 {x{t), P) = ^(^ + ^2,/3)-x(^,/3) ^ ^+(^(^)^ ^) ^ ^aA+z, 

^2 ^2 (1.3) 

VI {x{t), t, (3) = ^ b-{x{t), t) + -aA-z, 

Tl Tl 

which is properly defined because the diffusion shocks aA'^ z/t2 and aA~ z/ti are 
proper random variables with distribution /a+^ \'") ^ \ P \] ■ fact it is 



straightforward to show that 

E[v2 {x{t),t,P)\x{t)]^b+{x{t),t), 
E[vi {x{t),t,(})\x{t)]^b-{x{t),t), 

4 

Cov [vi\ x{t)] — Cov [U2I x{t)] = —acF 

T 



T 
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using the fact that E [t2^] = E [ti^] = 2/r from equation Hl.l|l . The expectations 
and covariances were calculated conditional on the collision occurring at time t. 

To define the solution to the difference equation properly, fix a time interval 
[0,T] and consider the collision set in this time interval x{tj,P),j = 1,...,N, for 
a random variable N such that < T and t^+i > T. Hence N is the last 
collision in the time interval [0,T]. Define N{t) — max^gjo i < t\ with 

tmn = tN(t),tmx = tN(t) + l then 

x(0), iiN{t)=oo, 

where = oo is an event that has probability zero as the particle inter-collision 
times are all finite and identically distributed. Hence P \N{t) < oo] = 1 so that 
tmx and tmn are well defined with probability 1. An alternative expression used 
by the Theorem below is the following martingale representation x(t,l3) = xn + 
b{xN,tN) {t~tN) + <jA{t — t]^) which differs from H1.4() only by the proportionality 
on the Gaussian shock. 

The point of the Theorem below is to show under what circumstances x{t, (3) — )■ 
a.e., where the position process x{t) is a strong solution to the stochastic dif- 
ferential equation corresponding to H1.2a|l with drift t) and variance matrix 
a — a{x{t), t). 

Theorem 1.1. Let the forward drift terms b~^{x,t) : M" x [0, T] i — > M" and vari- 
ance matrix a(x, t) : x [0, T] i — > M" x R" satisfy the following growth and global 
Lipschitz condition 

\\b+ {x, s) - b+ {y, s)\\ + \\<j {x, s) - a {y, s)\\ < K \\x - y\\ , 

\\b+{x,s)\\' + \\a{x,s)\\'<K'{l + \\xf), ^^-^^ 

for all x,y Cz M" and s G [0, T]. Let the drift function also satisfy a Lipschitz 
condition in time such that a constant Mf exists so that 

\\b+ {x,s)-b+ {x,t)\\ <Mf\s-t\, (1.6) 

for all X e R", and all s.t [0,T]. Let x{t,(3) be the position of the main particle 
at time t due to a finite set of collisions as specified in (|1.2a|l and (|1.3|l , i. e. 

x{t,(3) = a:o + ^ ^2 {x{tj,f3))Tj 

(1.7) 

= xo + Y^ [b+{x{tj,P), tj)Tj + crA+zj] , 

where r,- = tj+i - tj,0 < j < N{t), T^^t) = t - tN(t)-i, to = 0,tAr(4)+i = t (not a 
collision point) and where A'^Zj = z{tj^i) — z{tj), z{t) € R", j — 0, 1, N{t) are 
the Gaussian pulses. 

Let xo = x{(3, 0) = a;(0) be a random variable such that < oo then 

x{t,f3) x(t) almost everywhere as P oo for the process x(t) satisfying the 
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stochastic differential equation 
x{t) = Xq + 



b {x{s),s)ds+ / adz{t) 



(1.8) 



Remark 1.2. Condition (|1.5(l together with E[xq] < oo insure that there is a strong 
solution x{t),t £ 3?, to equation H1.8|l . This means that x{t) £ R" is a continuous 
process adapted to the filtration {g't; < i < oo} of the probability space (f^,^?, P) 
and that 



p 



[\b+ {x{s),s)\+al {x{s),s)] ds < 



oo 



1, 



(1.9) 



for all i,j = 1, ...,n. Condition (|l.t)|) is not required for a strong solution but is an 
essential condition for the convergence. The filtration {^t] < t < oo} encompasses 
the (T— algebra of events associated with the initial condition xq (a random variable) 
and the cr-algebra associated with the stochastic process z{t),t > 0. The solution 
is unique (any two solutions are equal with probability one) and can be written 
as a functional of the initial condition and a realization of the Gaussian process 
z(i);0 < t < oo. 

Remark 1.3. A strong solution for l|1.8|l implies a growth condition on the second 
moment of x{t),t > such that 

E[x{tr]<c[l + E[xl]')e^\ 

for some appropriate constant C > 0. In the case of few interactions no collisions 
occur before t = t^ (/? i 0) so then (|1.7(l reduces to 

x{t,0) =a;o+ [b+{xo,0)t + aA+zo] , 

with A+^o ^ z{t)- z{0). Then 

E [x{t, 0)2] = E [xl] + E [b+{xo, tjft^ + a^t] 

< E [xl] + K^E [xl] t^ + a^t < oo, 

from the growth condition (|1.5|l . 

Proof. For the sake of convenience the variance matrix will be assumed to be con- 
stant which docs not materially alter the proof. Fix N then 



\\xit,f3)-xit)\\ 



< 



j=N 

E 
E 

3=0 



tj+1 



[b+ {x[tj,(i),t^)-b+ (a;(s),s)] ds 



(1.10) 



|6+ {x{tj,(3),tj) - b+ {x{s),s)\\ ds. 



Because of (|1.5ll and (|l.t)|) the integrand can be majorized as 

{x{tj,P),tj)-b+ {x{s),s)\\ 
= \\b+ {x{tj,f3),tj) ~ b+ {x{s),tj) + b+ {x{s),tj) - b+ {x{s),s)\\ 
^ \\b+ {x{tj,(3),tj) - b+ {x(s),tj)\\ + \\b+ ix{s),tj) ~ b+ {x{s),s)\ 



< 



\b+ {x(t,,(3),tj) ~ b+ {x{s),t,)\\ + Mf \t, - s\ , 
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and applying this to Hl.l()|l results in 
\\x{t,^)-xit)\\ 



- X! / Wh"^ {x{t J, (3), tj) -h-^ {x{s),tj)\\ds 



f'^+' (1.11) 

^ ' -Ids ^ ' 



3=0 " *J 

3=N ^j. 3=N 

<KY, sup \\x{t,,p)~x{s)\\T,+^J2^3- 



j=0 *J<«<*J + i " j=o 



and so 



sup \\x{tN,(3) - x{s)\\ 

tN^S<.t 

j=N ^j. 3=N 

<kY, sup \\xit,,f3)~xis)\\T, + ^J2^3- 

j=0 J— — 3=0 

To extract the growth of the sup^.^^^f^^.^ \\x{tj , (3) — x{tj)\\ term from this in- 
equality the following discrete version of Gronwall's Inequality will be applied, see 
Shreve [TT]- 

Lemma 1.4. Let a„ > 0,7i = 0, 1, ... be a set of numbers such that 

{P,a,) + Qn, n=l,..., (1.12) 

3=0 



for positive numbers Pj,Qj, j = 0,1, .... Then for all n > 

U ( fc = "-l TT \ 

a„<aoPo^^+ 0«+ Y -^PkQA. n^\,..., (1.13) 
tin \ lli- / 



where n„ = 11^=0 (1 + i'j ) , « = 0, 1, .... 
Proof. From ((T~T^ it follows that 

3=n-\ 

a„ = ^ PjOj +Qn-Zn, n=l,..., (1-14) 

3=0 

for some positive sequence of numbers z„ > 0, n = 1, .... Let /3„ = X]j=o ^3^3 ^ ^ ~ 
0, 1, then p.l4|l can be written as 

/3„ = (1 + Pn) Pn-1 + Pn (Qn " ^n) , U = 1, 

with initial condition Pq — Pq^o- The solution to this equation equals 

n n 



Ho Hfe 
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using n„ = 11^=0 i^ + Pj)- Finally then 

PnCln — /^n Pn—1 

/TT -TT \ k=n~l 

A:— n — 1 

n„-i , ^ I ^ n„_i 



/ fe=n-l \ 



\ A;=l / 

since n„ - n„„i = n„_i (1 + Pn) - n„_i = P„n„_i. Now /3o = Pooo and 2:„ > 
0, n > 1, so dividing by P„ it follows that 

TT / k=n~l \ 

an < aoPo-^ + On + V ^^PfeQfe . (1.15) 

If all Pn,Qn, n = 0, 1, ... are equal so that P„ = P, (5„ — Q,n ^ 0, 1, ... then (|1.15|l 
reduces to 

an < (Pao + Q){1 + PT-\ n=l,.... 
Notice that the righthand side of p.l5|l increases monotonically with n. □ 

Applying itTT^ to |(lTT|| using n„ = 11^=0 (1 + = > and 

, , n "> n Tin^'n ^"rm lQc;^" r-om ci 

2 /1^J=0 'j 



Qn — Si-o — 0' with the last remark in the Lemma, it follows that 



sup \\x{tN, 13) - < ||x(0, 13) - x(0)|l Kto^ (1.16) 
„<.t<t tin 



j=N k=N j=k-l 

i=0 j=0 



Now a;(0, /?) = a;(0) and also 

j=N 

hence p.l6f) reduces to 

max sup \\x{s, P) — x{s)\\ 

0<j<N t.<s<tj + i 



j=N 



J=0 



SO that finally 



sup \\x{t,l3)-x{t)\\ 

0<t<T 



L 
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where 



L = ^MfT (1 + KTe^'^) . 



This is the result of the fact that E 
From Chebychev then 



Pr 



sup \\x{t,2''') - x{t)\\ > 4 



— 3Tt and the fact that ^Ar+i — t. 
L 



< 



2fc' 



so that by Borel-CanteUi x{t, (3) 
T I 0. This completes the proof. 



r;{t) almost everywhere on the interval [0,T] if 

□ 



Remark 1.5. This result justifies posing equation H1.2a|l to specify the discrete 
dynamics of the main particle and use the solution x{t),t > 0, (|1.8() as an ap- 
proximation of x{t,/3),t > 0. Clearly the discrete subordinate process that solves 
(|1.2a|l is difficult to determine. The theorem above does not address the validity 
of equation ljl.2b|l nor does it provide a hint as to the form of the backward drift 
b~{x,t). 

Remark 1.6. Unfortunately, the types of solutions that are most of interest have 
drift terms b~^{x, t) that can be singular and therefore do not satisfy conditions (|1.5f) 
and (|1.6() . However, Carlen showed in |2j that for a wide array of interesting drift 
functions a weak solution exists for x{t),t G M". Hence a filtration {^t; < t < oo} 
exists for which {z{t); < t < oo} is a Brownian motion and for which conditions 
(11.8(1 and (|1.9|l are satisfied. In this case uniqueness does not necessarily hold (in 
the same fashion) and no (measurable) functional exists to map the solution given 
the initial condition and the Brownian path. It is not clear at all how Theorem 
can be generalized to appropriate weak solutions of equations (|1.8|l and (|1.9(l . For 
more discussion on strong and weak solutions in the present context, see Shreve |llj . 
Rogers UHl or Carlen [121 . 

Following Nelson [Jj and Carlen ^2] it can be shown that the backward incre- 
ment for a continuous stochastic process is equivalent to a time reversed Markovian 
process with a related drift and Gaussian increment. Specifically the following 
applies. 

Theorem 1.7. If conditions (ll.5|l and apply, let a = al so that the variance 
matrix a is constant. Then, for a diffusion process x(t) G M" the forward motion 
of the particle x{t + T2) — x{t) and the backward step x{t) — x{t — ti) conditional on 
x(t) — X are given by 

x{t + Ta) - x{t) = A+2;(t) = b+{x{t),t)T2 + crA+z, (1.18a) 
x{t) ~ x{t - Ti) = A^x{t) = b^{x{t),t)Ti + crA^z, (1.18b) 

b-{x,t)^b+{x,t)-a^^P^, (1.18c) 
p{x,t) 

where p{x, t) € R,x € R" is the probability density for finding the main particle 
at position x at time t. The shocks A~^z and A~ z are independent Gaussian in- 
crements with mean zero and variances T2 and ti respectively. From this it follows 
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p{x,t)t = -V(6+(a;,i)p(x,t)) + A,p(x,i), (1.19a) 
pix,t)t = -V {b-{x,t)p{x,t)) - A,p{x,t), (1.19b) 



Pix, t), = -V ( + ^(^^ ^) , ^ (^^^9^) 



that 



so that 



which is the continuity equation. Here V = ^gf^j gf~^ ; = (^cf 
and (.)j denotes the time derivative. 

Proof. See references Nelson 7. and Carlen ^21 and the references therein for proofs 
where the cohision times are infinitesimally smalL For the justification of (|1.18a|l . 
(|1.18b|) and p.l8c|) see remark below. □ 

Remark 1.8. From this point onward no distinction will be made between x{t,(3) 
and x{t) as the latter can be obtained as the limit of the former. This is not only 
to relieve notation but also to emphasize the fact that the results that are to follow 
are only true in the limit of /3 t oo- 

Remark 1.9. Equation (jLlSafl is exact by definition and its solution approaches 
x{t),t > 0, which has a probability density p{x,t) and a forward drift &+(x,t). 
Therefore the backward drift b~{x,t) can be obtained as a result of (|1.18c|l . As a 
result equation (jl.l8b|) is approximately correct for small discrete time steps. 

The Corollary below uses Ito's Lemma and (|1.8|l to find an approximation for a 
function of the position x{t). 

Corollary 1.10. For any function 

A+/ = f{x{t + r2), t + T2) - f{x{t), t) 

= ift+a^AJ)T2 + \7f.A+x{t) + 0{T2) 

and 

A-/ - f{x{t), t) - f{x{t - n), t - n) 

= {ft - o^AJ)t^ + V/.A-x(0 + 0(Ti), ^^'^^^ 
where A+, A^,Ti and T2 are introduced in Theorem ()1.7fl above. 
Proof. Equation (|1.2U|) follows from Ito's Lemma, while for (|1.21|l it is clear that 
f{x{t - Ti),t - n) = f{x{t) - A-x, t - n) 

= f{x{t),t) - ftn - fxA-x + i/,, (A^x)' + 0(ri), 

so 

A-f = f{x{t),t)-f{x{t-n),t-n) 

= ftTi + f,A-x - i/,, {A-xf + 0(ri) 

= [ft - (T^fxx)ri + fxA-x{t) + 0(ti), 
which settles the proof. □ 
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Remark 1.11. Consider the one-dimensional case where G R and assume that 
the variance term depends on x{t) and t expHcitly so that a = (j{x{t)^t). If the 
variance process is at least once differentiable then the backward and forward rep- 
resentation can still be obtained as follows. Let y = f{x, t) = cr(p.t} '^^'^ 
Pf{y) be the density fmiction for y at time t then using Ito's formula 

A+y = F+{y, t)dt + ct,.A+z 



where 



ft + b+Vf ^ 
A^y ^ F^{y,t)dt + GrA" z, 

F-{y,t) = F+{y,t)^a'Pf^y^ 



dy 



and where A+z and A^ z are the increments as defined above. Inverting the func- 
tion above x{y) = f~^{y,t) and writing ^ — p^, the probability density for y can 
then be translated back. 

Returning to the backward and forward representation of the motion of the main 
particle in (|1.2a|) . H1.2a|l the following is now straightforward. 

Lemma 1.12. Let the position process x(t) € M" and the inter- collision times be 
2nd order Gamma distributed. Then at the collision time t 



E 



A+x{t) 



T2 



x{t) 



E 



E 



A+x{t) 



T2 



b+{x{t),t),E 
xit) 



A-~x{t) 



x{t) 



b-{x{t),t), 



x{t),tf+'^, 

T 



(1.22) 



A-x(t) 



x{t) 



= b-'{x{t),tf + 



2a' 



where t ^ 2/ (3 is the mean inter-particle collision time and e = a' jM is the 
diffusion per unit mass. 

Proof. The expectations are straightforward as the terms on the righthand side of 
the second and third equation in (|1.3|) arc finite curtesy of the fact that 



E 



A+zV' 


^ E 


' 1 ■ 


2 


T2 J 






T 




= E 


■ 1 ■ 


2 






.n. 


T 



E 

The equality in the first equation in l|1.21|l is due to the fact that 

' Px{x,t) 



ib+{x,t)~b-{x,t))pix,t) 



p(x,t) 



p{x,t) = 0, 



so that the average velocity of the main particle is the same whether a forward or 
backward view is developed. □ 
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Remark 1.13. As a result of (|1.22|l the forward and backward energy for the main 
particle can be defined as 

2 



A+x{t) 



T2 



M 



-E 



A-x{t) 



T2 



x{t) 
x{t) 



^b+{x{tltf+'-, 



— b-{x{t),t)^ + -, 



(1.23) 



because ^2a^/T = e/r. 



Remark 1.14. Recall that it was assumed that t ^ t otherwise ti is not sufficiently 
close to a 2nd order Gamma distribution. 

Remark 1.15. In the typical stochastic mechanics setting, — H/M so that the 
backward and forward energies have a drift component and depend on a con- 
stant equal to h/r. Typically in molecular applications at room temperatures 
Mb'^{x{t),t)^ » h/r and only at lower temperatures Mb'^{x{t),tY ~ h/T. How- 
ever in other applications like astrophysics or economics it remains to be seen what 
the drift energy is in proportion to the diffusion energy. In the molecular-kinetic 
theory of heat the diffusion coefficient equals cr^ ~ kT/rj-u so the mean versus dif- 
fusion energy ratio then depends on the temperature T and the viscosity 771, . 

The following example demonstrates a consequence of Theorem (|1.7() and equa- 
tions H1.2a|l and H1.18c|l . Consider a Brownian particle that arrived in x{t) — x 
and assume that the collisions occur closely in time so T2 and ri are small. From 
(I1.18c|) follows that the forward and drift for Brownian motion are determined as 
&+(x,t) = and b^{x,t) — —x{t)/t. The approximate energy gain or loss for the 
main particle can now be approximated from the outgoing and incoming energies 
around the collision point so in one dimension 



2{H+j{x{t),t)-Hl,{x{tht)) 



Xit + T2) -X{t) 



M 



x{t) — x{t — Tl) 



T2 J V Tl 

where {x{t),t) — H^,j {x{t),t) is the energy gain or loss due to the collision in 
x{t) at time t. 

Averaging over all collision points this reduces to 

2E[n+,{x{t),t)-nij{x{t),t)] 
2" 



ME 



a 

2M— 

T 



A+x{t) 



T2 



ME 



ME 



A-x{t) 



(1.24) 



2M— 



Ma' 



The expectation includes all paths emanating from the origin so this suggests that 
the main particle sheds energy continually. The expression does not depend on the 
mean inter-particle collision time or relates in any other way to the characteristics 
of the impacting particles. In fact the large contributions of the diffusion terms 
in the backward and forward momenta cancel exactly leaving the time dependent 
drift -MaV(2<). 
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Equation H1.24|l should have a more comphcated appearance to reflect the dis- 
creteness of equation H1.2a|l . However, for sufficiently small r the conclusion must be 
that a particle following a Brownian motion path must be flooding its surrounding 
heatbath with energy thereby raising its temperature. In the continuous represen- 
tation the energy transfer is inflnite at the origin which is an unrealistic physical 
anomaly due to the discreteness of equation (|1.2a|l . It is clear that a continual 
energy exchange is possible (viz. high energy protons fired into a plasma torus) 
but it is unexpected that a Brownian (main) particle is not in equilibrium with its 
surroundings. 

If not Brownian Motion then the question is what stochastic process describes 
the equilibrium motion of the main particle in a heatbath such that there is no 
energy exchange between main particle and heatbath? To answer this question 
the next Section shows that the main particle equilibrium motion is related to the 
total energy conservation and investigates the constraints on the probability density 
function for the position of the main particle. 



2. Energy Conservation 

This Section combines the properties of elastic collisions to show that the main 
particle plus heatbath particle kinetic energy equals a quadratic expression of the 
backward and forward velocities. From this it is shown that if the main particle 
is in equilibrium with the heatbath then there is only one acceptable probability 
density for the position of the main particle. Another conclusion in this Section is 
that the forward and backward velocities of the main particle through a collision 
point can only be independent if the heatbath particle motion through the collision 
is highly correlated. Examples are presented in the form of the Gaussian Wave 
packet and Brownian motion. 

To analyze the energy exchange associated with the collisions consider an elastic 
two-particle collision in n (typically 2 or 3) dimensions at time t between the main 
particle of mass M and the heatbath particle ni where typically {M > m). The main 
particle has pre- and post-collision velocities V2,vi £ M" while the heatbath particle 
has pre- and post-collision velocities W2,wi € R" respectively. By assumption the 
particles exchange energy and momentum during the collision hence V2 ^ vi and 
W2 wi. As there are no other interactions, the momentum and the energy during 
the collision must be conserved so that 

P2+ q2 = Mv2 + mw2 =pi+ qi^ Mvi + mwi, 

-(Mhl ) = __ + __ (2.1) 

= -(MK| +-KI ) = ^ + ^- 

To obtain a full solution to equation H2.1|l assume a linear relationship between the 
main pre- and post-collision particle velocities fi , f 2 G M" and the incident pre- 
and post-collision particle velocities Wi , if 2 G M" as follows 

where P = P(x, t), Q = Q{x, t),V = V{x, t) and G = G{x, t) are M" x M" matrices 
so that r = r (x, t) is a M^" x M^" matrix. The following theorem determines the 
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form of the matrices P,Q,V and G under the conservation of energy and momentum 
constraints H2.1|l . 



Theorem 2.1. Assume that the T matrix in (|2.2|) can be decomposed as T 
then the energy and momentum conservation in H2.1|l is equivalent to 

ml \ m 



P 

.V G ) 



(2.3) 



\ U ml \ U m 

or equivalently 

P^M + R^m = M, P^MP + V'^mV = M, (2.4a) 

Q'^M + S^m = m, Q^Af g + G'^mG = M, (2.4b) 

P^MQ + V^mG ^ 0. (2.4c) 

If M — MI,m — ml (I being the unit matrix) then these equations can be solved 
to yield 

U'^U = I, 

where U is an arbitrary n x n unitary matrix U^U — I . 

Proof. A straightforward but lengthy proof for this can be found in Appendix B. □ 

Remark 2.2. Notice that Q '^'^V, P + Q = I and V + G = / so the colhsion 
matrix r(a;, t) can be easily expressed in terms of the matrix Q. 

Remark 2.3. The matrix r{x,t) in H2.2|l consists of all center of mass collision 
information to translate the pre-coUision velocities vi , wi information into the post- 
collision V2, W2 configuration. In one dimension this matrix is naturally absent 
and in higher dimensions the matrix varies randomly from collision to collision. 
The matrix depends on F = F [x, t, Z) where Z is an asymmetric matrix Z = 
I — 2{I + U)~^ with U defined in 12.511 . The matrix Z will be referred to as 
the collision scattering matrix. With elastic collisions this matrix is independent 
of the energy pre-coUision energy vi and wi and typically reflects the physical 
circumstances of the collision event. 

For U = I, Z = 0, equation (|2.5|) reduces to a simpler set of linear equations 
which shall be referred to as the "simple" elastic collision 



V2\ (cos{9)I 7sin(0)/ 



; - 1 ; ' """Vflu I ( " ) = " ( ' (2-6) 

W2 J \ ^ I — COS(0)I I \Wi J \Wi J 

with obvious definition of the matrix Q and / being the unit matrix in n dimensions. 
Here 7 = cos(0) = , sin(6') = j^^- In one dimension this solution is unique 
as in that case Z = 0. This expression also applies if 7 = 1 (m = M) in which 
case the particles simply exchange velocities V2 = wi,W2 ~ wi • The matrix Q is 
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a reflection transformation with det(Q) 
eigenvectors 



cos I 
sin I 



with eigenvalues 1,-1. 



(2)/ \ ^^"V2, 
Notice also that cos(0/2) 



= -1, O^fi = /, Trace(ri) = and has 



(f) 



and sin(6'/2) = j^. 

Remark 2.4. Equation (|2.1|l has a unique solution in the sense that the main and 
heatbath particle exchange momentum and energy along the center of mass line 
at the moment of collision while other components of the motion remain invariant. 
This means that in higher dimensions equation (|2.6(l applies to the pre - and post- 
collision velocity components of W2,W2,fi,wi along the center of mass hne. The 
components of vi and wi orthogonal to the collision center of mass line remain 
invariant under an elastic collision. In this Section no real distinction is made 
between the properties of the heatbath particle and the properties of its center of 
mass line components. 

A simple Theorem now describes how the energy conservation for the motion 
for the heatbath particle m,wi (post-collision m,W2) and the main particle M^vi 
(post-collision M, ^2) can be expressed exclusively in terms of V2 and vi if 1)2. 6|l 
holds. 

Theorem 2.5. Let the momentum of the main particle and interacting particle be 
presented as pi = Mvi (post-collision p2 — MV2) and qi = Mwi (post- collision 
q2 = MW2) with vi,V2,Wi,W2 S M". Then the total energy TLt — ^(Mw| -\-mW2) — 
^{Mvf -\- mw\) is related to the pre 
particle as follows: 

1-1-7^ 

for the heatbath particle, while for the main particle 
M _ 



and post collision of momenta of the main 



1 


«2 + qi 


2 

H 


4 


q2 - qi 


2 


2 






2 



1 


P2 +Pl 


2 

H 


1 


P2-P1 


2 


2 






2 



W2 + Wi 



2 2 



1-f 7' 

In terms of the velocities these equations become 

74 ^ 1 
Mt 2 

Mt 2 

with Mt = M + m = M (l + 7^ 
momenta can be derived as follows 



W2 — Wi 



V2 -\- Vi 


' 1 


V2 - Vl 


2 


^ 272 


2 



(2.7a) 



Also a more direct relationship between the 



V2 + Vl 


' 1 


V2 - Vl 


2 


^7 


2 



= n l^ll 



\W2\ 



Proof. See Appendix A. 



(2.7b) 

□ 



Remark 2.6. The previous remark H2.4|l implies that in higher dimensions equations 
(|2.7a)l and (|2.7b|l apply to the center of mass line components of the main and 
heatbath particles. In other words if n > 1 the kinetic energy Hk incorporates only 
the parts of the motion of the heatbath particle that are altered due to the collision. 
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The motion perpendicular to the center of mass hne remains invariant hence this 
energy component remains the same and must be added to Tik in order to derive 
the full energy of the heatbath and main particle combined. In Section 4 equations 
(|2.7a|l and Ij2.7b|l will be investigated for the case where Z ^ Q so there TLk will 
represent the combined kinetic energy. 



Remark 2.7. It deserves mentioning that equations 12.7a|l and (|2.7b|l only hold for 
a collision point and are in fact incorrect for any intermediate points. 

Conceptually equation H2.7a|l suggests that the fracturing of the path of the 
main particle provides a direct indication of the amount of energy that is involved 



in the main and heatbath particle system combined. Moreover, condition (|2.7b|l is 
independent of the energy constraint H2.7a(l and tends to constrain the momentum 
exchange between the particles involved. For this reason H2.7a|l is referred to as the 
"energy" balance/constraint and (|2.7bl) as the "momentum" constraint. 

The agenda for further developments is now to employ the collision representa- 
tion of the Markovian process describing the path of the main particle in equation 
l|1.2a|l and substitute the respective forward and backward velocities (11.3(1 into 
(|2.7a|l . Due to the distribution of the inter-collision stopping times the expectation 
of the total energy components in the collision process will be well defined. The 
mean inter-particle collision is small so the use of equation l|1.2b|) is well justified. 
The total energy as a function of time and position can then be investigated. 

First an important consequence of equation 1(2.2(1 must be noted. Rewrite equa- 
tion 1(2.6(1 so as to represent the pre- and post-collision velocities of the main particle 
as a function of the velocities of the colliding heatbath particle. This shows that 



V2\ _ 7 (cos{e)I I 



vij sm{6) \ I cos{9)I J \wi 



W2 



(2.8) 



hence the statistical properties of V2,vi are generated by the behavior of W2,'Wi 
and vice versa. 

Specifically equations l(1.2a(l and 1(1. 2bp demand that the motion of the heatbath 
particle must be of the following form 

A+ 

T2 
Tl 

which substituted into equation 1(2.8(1 yields 

A+z\ / A^z 
b'^{x{t),t) +a =cos(6') b~{x{t),t) +<T 

T2 J \ Tl 

+ 7sin(^) (^g-ixit),t)+Lo^ 

+ / . X ^ A+\ sin(6l) /, , , , , A'z 
g+{x{t),t) = ^ b-{x{t),t)+a 

T2 / 7 V 

-cos{9) L'{x{t),t)+uj^ 
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or equivalently 



b-{x{t),t)^ 



7 /cos(6l) 1 ' 
sin(6i) 1 cos(e')^ 



7 / cos(0) 
sin(6') V 1 



g-{x{t),t) 



cos{0)) L^, 



(2.9a) 
(2.9b) 



Hence the motion of the hcatbath particle must be driven by backward and for- 
ward drifts g^{x, t) and corresponding random impulses lu^^, uj^^. As subsequent 
collisions involve different heatbath particles the realizations from the hcatbath par- 
ticles wi, 1^2 should be the result of a time dependent random field. Specifically 
the post-collision heatbath velocity W2{x{t),t) for the collision at x{t),t is not equal 
to the pre-collision velocity wi{x{t + T2), t -f T2) as they involve different particles. 
No continuity equation therefore applies to g^{x,t). However equations H2.9a|l and 
(|2.9b|) show that the heatbath realizations are constrained by the motion of the 
main particle. 

The other important consequence of (|2.9bp is that the main particle forward 
and backward velocities A+z and A~z can only be independent if the heatbath 
particle forward and backward velocities A+ and A~ are appropriately correlated. 
The following Theorem details the autocorrelation in the heatbath particle motion 
and also shows that the random terms A+ and A^ must be cross correlated with 
A+z and A~z. 



Theorem 2.8 

in equation 



If the main particle follows a Markovian path i.e if A'^ z and A z 
~ '2 - T^^, where = 7^/(1 -I- 7^^). Also 



I2.9bp are uncorrelated, then uj 
the A+ and A~ must be correlated as 



corr 



T2 



1 - 2a' 



If m << M then a is very small hence the variance term lo must be large and 
the correlation between the pre-collision and post- collision velocities of the heatbath 
particle must be very high. 

In addition the forward and backward velocities of the main and incident particle 
have covariances 



E 



E 



Proof. Inverting equation (j2.9b|l yields 



[AiA+zl 




[A-A-zl 


[ r| J 


= E 


[ rl \ 








[A+A-zl 


= E 


[A-A+zl 


T2T1 




r2T\ 



V2c 



7sin(6') 
V2a 
7sin(0) 



cos{9), 



7sin(^?) 



- cos(6l) 
1 




(2.10) 
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and if the increments a^—^ and cr-^— ^ are independent then 















m 













2a' 



T72 sm\e) 



cos{e) 
1 



(7 



1 



1 - 2a' 



1 

- cos(e') 

- (1 - 2a') 
1 



where a' — j'^ /{I + j*). If it is assumed that 



E 



T2 



E 



Tl 



then Lo = a I (aV2) and the correlation foUows directly. 
Using H2.1()(l again it is also clear that 



LuaE 




T7 sin(6') 



so that 



E 



' /a+A+z 


A + A-z 




T2T1 


1 A-A + z 


A-A-z 






V2a / 


- cos(6i) 


jsm{9) [ 


1 



1 

cos(6l) 



which concludes the proof. 



□ 



Remark 2.9. If 7 << 1 then the uj' variance term becomes very large and the 

A^ A^ 

correlation between <J-f^ and c-^q^ becomes tighter. In the limit where 7 ~ 0, e.g. 
a football interacting with molecules or a planet interacting with cosmic particles 
(or photons) a small diffusion constant a for the main object is associated with an 
enormous heatbath particle momentum proportional to a/ (aV2). In this case the 
post-collision heatbath particle velocity size is identical to the pre-coUision heatbath 
velocity but moving in the opposite direction. 

The next issue to investigate is the total kinetic energy Hk in H2.7a|l as a function 
of time while the main particle travels from collision to collision. By assumption 
the collision is elastic and conserves the total energy so if wi and 'W2 refer to the 
same particle (the first one before the collision and the second one after the collision 
with the main particle) then by definition 



nk = nk{xit),t) =. ^M\v2\^ 



1 I |2 

-m |W2| 



-m\wi\ . 
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However, for the subsequent collision the main particle moving with momentum 
Mv2 particle will meet a different heatbath particle at time t + T2 in a new position 
x{t + T2) = x{t) + b'^{x, t)T2 + crA+z(t)/T2. This new heatbath particle will have a 
new momentum mq' which is chosen randomly from the ensemble. Obviously the 
exiting total kinetic energy at x{t) and the combined "new" total kinetic energy at 
x{t + T2) will not be equal since 



H'^{x{t+T2),t + T2) = \M\V2\' 



Im\w2\\ 



Hence that the combined system (main + heatbath particle) gains or looses 
energy ATi.k{x{t),t) equal to 

AHk{x{t),t) = Hkixit + T2),t + T2) - Hkixit), t) 
On average the amount of energy exchanged equals 



E[Ank{x{t),t)] 



m 



E 



1/2 



-E 



\W2\ 



(2.11) 



where the expectation ^E 



is the average energy of the new incoming heat- 



'-E 



\W2\ 



is the post collision energy 



bath particle at collision time t -\- T2 and ^ 
of the previous heatbath particle at the previous collision time t. The expectation 
here is over all paths, heatbath interactions and all positions x{t). 

If 1)2. 11|) is positive the original post collision energy of the heatbath particle at 
t has a lower energy than the new heatbath particle colliding at t + T2 and if this 
quantity is negative the collision accelerated the heatbath particle coming out of 
x{t) in comparison to the typical heatbath particle. In the first case the colliding 
heatbath particle returned to the heatbath with less energy than the average heat- 
bath particle and in the second case the collision accelerated the colliding heatbath 
particle beyond the heatbath average. If (|2.11|) is positive the heatbath puts en- 
ergy into the main particle lowering its temperature and if (|2.11() is negative the 
main particle radiates energy into the heatbath thereby heating it up. Changes in 
the expected total energy along the path of the main particle therefore show the 
amount of energy that is being exchanged between main particle and heatbath. 

The important result is that if the main particle resides in heatbath equilibrium 
the expectation H2.11|l should be zero. In other words if the main particle does 
not radiate energy into the heatbath, if the main particle movement has adopted 
the heatbath temperature then the total kinetic energy Tik must be a conserved 
quantity. If radiation occurs in the form of an energy exchange between main 
particle and heatbath then there must a potential or explanatory term describing 
the additional physical process. In this case it must be assumed that the total 
kinetic energy and the potential term together are conserved. 

Assume therefore a static potential <I>p G K, $p : R x [0, 00] 
(average) energy exchange can be related to (|2.11l) as 

E [<^p{x{t + ra), t + T2)] - E [<^p{x{t),t)] 



such that the 



2 



\W2\ 



(2.12) 
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then 



E [Hkixit + T2),t + T2) - Hkixit), t)] 

= -{E [%{x{t + T2),t + T2)]- E[%{x{t),t)]). 

Defining the total energy of main and hcatbath particle as 'HT{x{t),t) = Hk{x{t), t) + 
^p{x{t),t) then implies that E [Hrixit + T2),t + T2)] = E [nT{x{t),t)]. In 

other words the Hamiltonian E [Hxixit), t)] must be time invariant. The expecta- 
tion E [Hkixit), t)] is the kinetic energy term associated with the motion (of both 
the main and incident particle) and $p is the static potential energy term that 
regulates the energy exchange between heatbath and main particle. 

Using Theorem (I2.5|l , definition 1)1. 2a|) and approximation l|1.2b|) it is possible to 
provide more detail on the precise form of the total energy total energy Tix- 

Proposition 2.10. Let the average velocity of the main particle be defined as 

b+ixit),t) +b-ixit),t) 



vix,t) 



(2.13) 



then in n dimensions the expectation of the total energy E [Ht] can be written as 
E [Ht] --E[nk + $p] 



Mt 



-E 



\vixit),t)\' 



E[%ixit),t)] 



2ne 
sin^ (6')r' 



E 



(2.14) 



with Mt — M + m, — jj. Here pix,t) is the probability density function for 
xit) and the potential $p is defined in H2.12|l . 

Proof. The terms that requires an explanation are the expectation of the kinetic 
energy term Tik and the resident constant. Using H2.7a|l this expectation reduces 
to 

2' 



E 



-E 



1 



-E 



V2 + Vi 



V2 - Vi 



vixit),t) 



1 A+z 
—a 

2 T2 



1 



1 A^z 
—a 

2 ri 



1 



:E 



272 

which can be simplified to 



b+ixit),t) - b-ixit),t) 



1 A+z 



1 A-z 



E 



Ti-k 



Mt 
1 



^Ie 
2 



\vixit),t)\' 



:E 



b+ixit),t) - b-ixit),t) 



= '-E 
2 



vixit),t)y 

n (1+7^) e 
2mr ' 
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^ 1 + ^ — 



2 / ^ 
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(2.15) 
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if cr^ = jj. Again here p{x,t) is the probabihty density for x(t) used in equation 
()2.7a|l to find an expression for the &+ ix{t), t) — b~ {x{t),t) term. 
The constant in E [Hk] /Mt is the resuh of the fact that 



(l + 7')e _ (1 + 7 



a constant depending only on the properties of the incident particles. This means 
that 

E [Ht] ~ nMT - 27^7 ^ 

which explains the constant in H2.14|l . This reconciles all the terms and the Theorem 
is proved. □ 

Remark 2.11. Comparing the diffusion energy term e/r in (|1.23|) it is clear that the 
diffusion energy for the total kinetic energy ^^.^j^g^a^ in Proposition H2.10|l above is 
much larger if 7 << 1. This is due to the presence of the heatbath particle diffusion 
energy. 

The radiation requirement introduced above imposes a restriction on the distri- 
bution density for the position of the main particle. This result is summarized in 
the following Theorem. 

Theorem 2.12. Let the potential $p : R x [0,oo] i — * R defined in equation (I2.14|) 
be such that 



d_ 

It 



-E[% {x{t),t)]^E 



b+ + b- 



.V0 



(2.16) 



for a suitable potential : K. x [0, oo) i — > M which is at least once differen- 
tiable. Then the only probability density distribution p{x, t) for the main parti- 
cle position process x{t) that allows the total energy E [Ht] — E [Hk + <&p] de- 
fined in Proposition (|2.1U|I to be time invariant is derived from the wave function 
ip{x, t), a; G R, t € [0, oo), ■0 : R x [0, oo) i — > C such that p{x, t) = |V'(a;, t) where 
the wave function satisfies Schrodinger's equation 

ixi^ix, t)t = A,^(x, t) + (t>{x, t)^{x, t), (2.17) 

with X = MtV = Mt^ ^ i) e = 2e/sin(0) and A, = // 

ri{x ,t)+iS{^ ,t) 

the wave function is written as = ^p{x,t) — e x then the forward and 

backward drift can be written as 

b^{x,t) = ^{VS±jVR) 

= i^(Im±7Re)MM ^''''^ 
The constant total energy now equals 

E [Wt] - [l VV'fl + E [%] + (2.19) 



2Mt 



sin^ (6l)i 
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Proof. The proof here is part of a more general result presented in Appendix [HI 
Nelson [H] showed the relationship between the wave function and an energy func- 
tional similar to H2.14|l . In Nelson T this result was extended to incorporate a 
potential term as introduced in H2.12() though condition (|2.16() is different. The 
proof in Appendix El is for a more general energy expression incorporating the 
presence of the collision scattering matrix Z but runs similar to the presentation 
in Nelson (Zj. Carlen ^21 demonstrated that the stochastic differential equation 
H1.8fl admits a weak solution if the potential <j>{x, t) belongs to a class Kato-Rellich 
potential. □ 

Remark 2.13. Equation (|2.19|l above incorporates the total energy of the system, 
i.e. the energy of the main particle as well as the energy of the incident particle. 
The mass in the equation refers to the combined mass of the system Mt = M + m 
rather than the main particle and it is not specified how this energy is distributed 
between the two particles. 

Remark 2.14. Under conditions H1.5|l . H1.6(l and sufhciently small t the discrete 
collision process (|1.2a|) can be approximated by the continuous stochastic differen- 
tial equation H1.8|l . Hence both the momenta and (forward and backward) energies 
(|1.3() can be approximated using (|2.18l) see Carlen 13^. As a result Ht in (|2.14|l is 
now a well defined estimate of the combined energy. 

There is no reference in this approach to quantum mechanics or a stochastic 
interpretation of quantum mechanics. Though the proof in Appendix IbI is quite 
similar to the stochastic mechanics approach presented in Nelson |7], Theorem 
(|2.12|l above does not reproduce an interpretation of quantum mechanics at least 
not perfectly as the scaling here is different. Obviously, if Planck's constant were 
chosen as h = x then (|2.17|l reduces to the wave-function for a particle, however, 
the forward and backward drifts in H2.18|l explode if the mass ratio 7 were allowed 
to approach zero. In addition the mass term Mt incorporates information on both 
the main and heatbath particle and the energy H2.19|l refers to the combined kinetic 
and potential energy. 

2.1. The Brownian Motion. This continued example demonstrates what hap- 
pens to the total kinetic energy for a particle following a Brownian motion if no 
conservation law applies. The Gaussian distribution is the result of the diffusion 
equation with a zero forward drift b'^{x, t) so that the probability density and back- 
ward drift equal 

p{x,t) = 

then equation ()2.15f) reduces to 




1 + 7^ \ Mtc 



., , 1-f 7^ \ Mre 
' 872^ ) 2rm 
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SO the total energy in classical diffusion decreases continuously though there is no 
evidence of a potential. Applying H2.11|l shows that 

E[Ankix{t),t)] 

= E [Hkixit + T2), t + T2) - Hk{x{t),t)] 

1 r 



7 sin(0)n 



E 



{t + r) t 



7sin(6')nt^ ' 

which implies that a particle following a Brownian path always radiates energy 
into the heatbath. Brownian motion seems a non-equilibrium solution and a main 
particle can only perform a Brownian path due to a hidden potential that forces the 
motion (a dampened time-dependent oscillator would suffice). The effect mitigates 
quickly with time t increasing as a multiple of the collision time. In a very dense 
heatbath the inter-collision time shortens so for macroscopic objects the moment 
to equilibrium must be almost instantaneous. Notice that the effect depends on the 
mass ratio and is not just a function of the main particle mass M . 

To understand why for a Brownian Motion the average total energy changes 
consider that classical diffusion is a limiting case (limit in time and space dimension) 
of a particle stepping forward over a lattice grid with equal probability. If the 
process has moved a large distance to position a; in a short amount of time the 
stochastic path of the particle has a strong backward drift in the direction of the 
origin. However, the forward drift is zero so for these paths the contribution to the 
total energy is large and ultimately too large to limit the average total energy. A 
path produced by a time invariant average total energy realizes that the backward 
drift is large and adapts the forward drift to point in a similar direction as the 
backward drift. Changing the drift will then reduce energy in the regions where x 
becomes large by reducing the acceleration. 

For the next example the following Corollary will be useful. 

Corollary 2.15. The drift for the heatbath motion is given by 
,±(.,t) = ^(v5T^Vi? 

^ ^ImT^Re^^^^^'*) 



Mt\ ^7 ) i^{x,t) ' 

Proof. Using equation H2.9a|l the relationship between the heatbath particles and 
main particle can be written as 

'g+{x,t)\ _ 1 ^~cos(g) 1 \fb+{x,t)' 



so clearly 



g-{x,t)J 7sin(6l) V 1 - cos(6l)y \b-{x,t) 

1 f- cos{e) 1 \ fVS + 7Vi?\ 
' MT7sin(6') V 1 -cosie)) [yS-jVRJ ' 

+ 1 ({1- cos(g)) (V5) \ 

^ ^ jMTsmie)\-l{l + cosi9))VR) 

Mt\ 7 J 



26 



JOHAN G.B. BEUMEE 



and the Corollary is proved. 



□ 



Remark 2.16. Analogous to Lemma (I1.12f) and equation H1.23|l the average momen- 
tum and energy for the heatbath particle can be determined as 



E 



T2 



X{t) 



'Hl{x{t),t) ^ -E 



g+{x{t),t),E 



A+wjt) 

T2 



A-w{t) 



Tl 



xit) 



9 {x{t),t), 



ix{t),t) 



-E 



A-w{t) 



Tl 



xit) 
x{t) 



2 Qf^r 



(2.20) 



-g-{x{t),t)^ 



7 e 



since mcr^/(2a^T) = 7^e/(a^r) 



This Corollary together with 12. 8|) completes a description of the behavior of the 
main and heatbath particles. If the main particle motion behaves like a martin- 
gale where its energy exchange with the heatbath is derived from a potential <f>p 
then its probability density and (forward and backward) drifts are derived from 
Schrodinger's equation (|2.17|l and equations H2.18|l . The colliding heatbath particle 
has a drift as specified in Corollary H2.15|l and exhibits a high degree of correlation 
between its backward and forward momentum as shown in Theorem (|2.8|l . The 
heatbath particle motion is also highly correlated with the motion of the main 
particle and can not be a martingale process itself. 

The following example shows a solution to equations H2.17|l and H2.18|l for the 
case where the path of the main particle can be represented by a generic Gaussian 
process. 



2.2. The QM Wave Packet. Applying Theorem and Corollary lfTT3|l to 

the Gaussian function wave packet in one dimension shows that the wave function 
can be represented as a Gaussian superposition of single momentum solutions to 
the wave equation H2.17|) . Hence 



OO 

I e ""^e 

-oo 



JztJ.{x,t)-^r{t) 



where 



1 t 



xMt 
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Hence 



where 



Zr{t) = T{t)T{t) = + "'i') , « = (2.21) 



-1 ,- Pn' 



and where C" = ( yT(t)xf'V\7^ ) c'^^mt ^ xhis constant is chosen to insure that 



p{x,t) = |^(a;,f)|^ is a probabihty density. 



I arccos 



Rewriting T{t) = y/Z^e K-IV^^J the wave function ^p{x,t) can be 

-R.(:c,t) + iS(x,f) 

represented as = ^{x,t) = e ^ where 

--Bf)' 1 



at ( Pot\^ , f XPQ Pot 



1 / 1 

- arccos 



2 \a^^/Z^)J' 
with a and Zr{t) as defined in H2.2f|l . FinaUy, then 



— i?(x,t) 



dx ' ' ^ x^GiZy{ty 

d , . at ( Pot\ ^ 

Hence x{t) is a (Gaussian) random variable such that £'[x(t)] = ^^,var{x{t)) = 

Zr{t)x'^l 

Using the drift funct ions become 

= ^rr^TiT # V"^ T (2.22) 



xAfTZr(t) V Mt J \ (JeJ Mt 
and the total energy becomes 



2 



, Po , 2e 



(2.23) 



2AfT 2Mt sin(6')2r' 

Hence the energy consists of a contribution due to the energy dispersion controlled 
by (Te, the mean kinetic energy term with po a-nd the diffusion term proportional 
to the variance ct^ 

Now the forward and backward energy ((T^ equal 

+ , ^^ Mai ( iV ' 



(.(t), i)] = ^^^^p^ [at ± -|J + + (2.24) 
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SO that there is a negative energy transfer E [h'j^j {x{t),t)] — E [HJ.j {x{t),t)] w 
~T^^^^f(t) 1^ ^ which decreases with time. This shows that the heatbath is ab- 
sorbing energy generated by the main particle initiaUy but this approaches zero 
depending on the size of the term a = 1/{x^t)- For large objects 7 w so then 
this heat loss for the main particle will be negligible. Since the total energy is 
conserved this argument suggests that the heatbath must be gaining energy as a 
function of time which indeed will be shown below. This is an example of a case 
where the total energy is conserved as demanded by the non-radiation condition 
but only because both component energies change in time. To conserve E [Ht] is 
therefore not equivalent to demanding that the backward and forward energies for 
the main particle are equal. 

The average energy for the main particle equals 

E[nM,avg ix{t),t)] = ^E [n+ , t) + 7^^, ixit),t)] 

Mai / 9 9 pI e 



In the case that at « - equivalent to a very short time step t « or a very small 
mass ratio 7 « - the average main particle energy reduces to 

E[nM,avg ixiO),0)] = ^E ix{0) , 0) + Hm {x{0) , 0)] 

M72 n2 g 

2M2 ^ r 

On the other hand from 12.2411 it is clear that 

E[TiM,avg (a;(oo),oo)] = ]^E [7Y|^(a;(oo), oo) + 7Y^^(x(oo), oo)] 

Mai Pi ^ ^'-''"^ 

= + M— h - 

2M| ^ 2M2 ^ r' 

1. Notice that the terms Mal/{2M^) and A/p^/(2M|) in the limit 



energies E [HM.avg {x{oo), cx))] and E [T-lM,avg (a;(0), 0)] are close in size to the first 
two terms in (12.23(1 as long as 7 is small. However the diffusion term is (|2.23|l is 
very large in comparison to the diffusion term e/r in E [T-lM,avg {x{oo),oo)] and 
E \HM,avg (a;(0), 0)] see the discussion below. 

To investigate the heatbath apply Corollary (|2.15l) to determine that for the 
heatbath particle 

1 ^ MU^,^ M,z^L 



so that 



" " xMrZrit) V Mt ) ^ lal J ' Mr' 



r M . . . .1 mat ( 1 \ Pn vauj 
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and 



E[H,n.avg{x{t),t)] ^-E [n+{x{t),t)+n+{x{t),t)] 



ma^ / 2 2 1 \ 



,2 



2A/2 



+ 



T 



Hence 



E [H„{x{Q), 0)] = -i? [n+{x{0), 0) + 7^+ (a;(0), 0)] 



(2.27a) 



Mai Po "^^^ 




and for large time t the average heatbath energy becomes 



E ['Hm,avg{x{oo), oo)] = -E (a;(oo) , oo) + 7Y+ (a;(oo), oo)] 



(2.27b) 




which indeed decreases from Mal/{2M^) to mal/{2M^) as time progresses. 

It is interesting to compare the variance contributions to the energies (|2.23|l . 
fll^ and ifT^ . Recall that mw^/r = mcr^/ (2a2). If 7 << 1 then mLj^/r >> 
e/r so the diffusion energy contribution to the heatbath particle is much larger than 
the diffusion energy for the main particle. The energy due to the dispersion term 
is at first a small contribution in the main particle energy in H2.24|l but then increases 
(see (|2.25all . (j2.25b|l ) while this term in the heatbath particle energy ()2.26|l does 
exactly the opposite (see (|2.27a|l . (I2.27b|) 'l. Both the main and heatbath particle 
energy are proportional to the same kinetic energy contribution weighted with their 
respective masses. This implies that the kinetic energy is almost exclusively carried 
by the main particle and for small 7 very little kinetic energy is carried by the 
heatbath particle. In fact in the small 7 limit the heatbath particles are all moving 
through the heatbath with energy Hm.avg ~ cr^/ (2q:^t) « e/ (27^t) without being 
affected by the presence of the main particle. After the collision the heatbath 
particle emerges with exactly the same speed and opposite direction. 

A curious consequence of this example seems to be that the statistical charac- 
teristics of the heatbath particles are affected by the energy constraint in a similar 
fashion as the main particle. The time dependent term in (|2.26f) shows a behav- 
ior that reflects the main particle and H2.26|l carries a kinetic contribution that is 
is also present in (|2.24l) . This is slightly unrealistic as it seems to imply that the 
heatbath particle behavior depends additionally on the main particle kinetic energy 
rather than external factors alone. Though the effect becomes very small as the 
mass ratio 7 decreases the only way to render equation l|2.26|l time-independence 
is to assume that the average collision time r changes as a function of energy or 
correlate the drift dependence between pre- and post collision velocities. The next 
Section addresses these constructions. 



The example in the previous Section showed that a constant total kinetic energy 
can be achieved but the individual main and heatbath particles energies display time 
dependent behavior. This Section attempts to investigate the momentum constraint 



3. The Minkowski Invariant 
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(|2.7b|l to impose energy conservation embedded in the heatbath particles through 
the point of coUision. It will be shown that the relationship between the inter- 
particle collision time and the total energy leads to a type of geometric invariant. 

To obtain an equation similar to H2.14|) for the momentum constraint substitute 
l|r^ into (|2Jb|) to find 



+ — jE 



2 



v{x{t),t)\^ 
b+{x{t),t) - b-{x{t),t) 



(3.1) 



2a2A/T' 



where 5+ = b^{x{t),t) and = b^{x{t),t) are the usual forward and backward 
drifts respectively, v{x{t),t) = {b'^{x{t),t) + b^ {x{t),t)) /2 as defined in (|2.13(l and 
a — 7**/ (1+7"*). Assuming that the underlying process is a martingale and applying 
equation H1.18c|l reduces this to 

l„r, .9 . ,,1 1 
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(3.2) 



The constant in this expression is the result of the fact that in one dimension 
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and multiplying with n for the multi-dimensional case yields the value for the 
constant in 1)3.2(1 . 

In a heatbath where the main particle is in equilibrium with the incident particles 
it should be expected that the post-collision energy for the heatbath particle is equal 
to the pre-collision when averaged over all paths and positions. Otherwise there 
will be an average heat loss or gain for the main particle. So if the main particle 
is in equilibrium it is expected that E [w^] = E [wf^ = where c is the average 
velocity of the heatbath particles. This means that i (i? [uj^] + E [wf]) — and 
defining A 



\AXrf = 



E \v{x{t),t)\' 



^E 



b+{x(t),t)~b-(x(t),t) 



(3.3) 



reduces equation (|3.2|1 to 



Hence 



2a2A-f7 



2—2 
C T 



2a2M' 

The results can now be summarized in the following theorem 



(3.4) 
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Theorem 3.1. Let x{t) G M" he the coordinate process of the main particle and 
define |Axt| as in equation Assume that the main particle is not radiating 

energy or receiving energy so that the backward and forward velocities of the heatbath 
particle are equal, e.g. E [uj^ — E \_wf\ = (? . Then the drifts must be such that 



,2=2 



C T 



\Axy-\ forms an invariant of the motion, i.e. 



3=N- ^ 
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ncT 
2a^M 



= constant, 



where is the average number of collisions such that E X]j=]^^ 



(3.5) 



Nrf = T. 



Proof. The proof is simply that (|3.4|l imphes that 

3=N- , j=N- 



El 
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ner neT 
2a^M ~ 2a2M' 



3=0 j=Q 

since is the (average) number of inter-particle collisions between time and 
time T. The righthand side is independent of the average inter-particle collision 
time. □ 

2 



Remark 3.2. If |Aa;-(, 



and is large then the solution to equation (|3.5() 



equals 



so there is a reference average inter-particle collision time for a slow-moving (sta- 
tionary) main particle in the heatbath. Combining this and H3.5|) shows immediately 
that 
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2=2 
C T 
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c^Tto. 



To gain some insight into this Minkowski type relativistic invariant H3.5|l consider 
the example of the Gaussian wave packet (|2.2|) in the previous Section. From 
equations H2.22|l it is clear that 

Pot 



1 at 
-(6+(.,t)+6-(.,t))=-^^^-^ 
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- (b+(x,t) - b~(x,t)) = - , , ^ 



_ Pat \ 
Mt) 



Po 
Mt 

0, 



Po 
Mt'' 



since at/Zy^t) « and 1/Zr{t) « for sufficiently large time t. As a result equation 
113.511 becomes 




(mt) 



eT 
a^M' 



This suggests that the average inter particle collision time r depends on the mean 
motion of the main particle specifically r ~ (c^ — {po/M)"^^ . However this pre- 
sumes that the drift \Ax-\ does not depend on the average inter-particle collision 
time T = t(|Axt|) which is unlikely to be reasonable. In general the mean inter- 
collision time and the drift will depend on each other and the most straightforward 
approach is to explore a linear relationship. 
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There does not seem to be a separate mechanism for introducing a dependent 
inter-cohision time such that H3.5|l holds except possibly a statistical correlation 
between the forward and backward velocities IS.'^ z and /Sr z. This creates a separate 
relationship between the correlation corr = Pa+z,a-zI {I being the 

unit matrix) of the coordinate process x{t) and the inter-particle collision time r. 
Returning to H3.1|l and introducing the correlation term shows that 



-E 



\v{x{t),t)Y 



E 



1 



\b+{x{t),t)-b-{x{t),t)\ 



A+z A^z 



IAxt 



IAxt 



T2 

o^Mt 
ne 



874 



vE 



A-z 



T2 



(3.7) 



2Mt 



( 9 \ 



if p\ = {\ — 2a^) Pa+zA-z- Multiplying with and dividing by (l — p.y) this can 
be written as 



2 2 
c r„ 



lA^J 



ncT 



where 



|A^^ 



£; \v{x{t),t)\' 



^E 



This means that (|3.5|l reduces to 

Y: {c'rl-\Ax..'' 



neT 
2a^M ' 



(3.8) 



where now the summation ranges over rather than over N^^. To obtain an 
estimate of the size of the correlation assume that r « to, write = \Axr^ |^ /t'^, 
then substitute H3.6|l into (|3.8(l to yield 



2—2 

c r,, 



E ( 

which is equivalent to 



I Ax J 



T 

To 



■ 2—2 

c r,, 



2— 2\ 



neT 
2a^M 



C TqT, 



1-4' 



(3.9) 



As a result then equation H3.7|l shows immediately that « jt? . 

Interpreting r and | A differentials equation 1)3.8(1 becomes the well known 

Minkowski invariant in Relativity Theory and its solution is the (linear) Lorentz 
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transformation which relates elapsed time versus motion as perceived in different 
reference frames. The following Theorem smximarizes the results and shows the 
simple linear transformation. 

Theorem 3.3. Let the main particle in the heathath he in energetic equilibrium 
with the heathath so that E = E \_wf\ = at any collision point t. Now let 
the forward and backward velocities A+z and A~z he statistically correlated so that 

Then condition 



corr 



^) - PA+zA-zI '^rid let pI = {I- 2^2) p^^^^_ 



1)3.1(1 is equivalent to 



3=0 



neT 
2a^M ' 



(3.10) 



where N— ~ T/t is the average number of collisions in time T assuming inter- 
collision time T and where 

T 



7a' 



\Axv\ 



E \v{x{t),t)Y 



^E 



b+{x{t),t)-b^(x(t),t) 



Now if V — Axu/r^ and v' = Ax/t the solution to equation H3.10|l equals 



Axv 



1 



\v — v' I" 



^2 



Axr 



(3.11) 



Proof. For the more formal solution to H3.10|l let Ax- be the drift associated with 
the inter-collision time r and Ax^ be the (larger) drift associated with the inter- 
collision time Ty then the invariance condition l|3.10|l suggests that 



Axv 



A B 
F E 



Ax-, 

T 



(3.12) 



for a constant matrix A G M"^", B. F'^ £ 
this into ((TTUIl yields 



_ / AAxr + Bt^ 
~ \FAx- + Et^ 

"^^ and E a constant. Substituting 



2—2 



c^ [Et + FAx-Y ~ \AAx- + Bt\ 



= (c2^2_ |^|2 



- Ax^ {A^A - c^F^F) A 



(3.13) 



where c^EF^ = AB is chosen to avoid mixing terms. Applying r = tq then 
Axto = so that (|3.13|) reduces to 



Ax^ 



and define B/£; = V = A^„/t„ e M". 

From c^EF"^ = AB it then follows that = Av/c^. Equation (pDUl now 
reduces to 



2—2 



Axt ^ E' 



2—2 
C T 



(3.14) 



- Ax^ {A^A - c^v'^F'^Fv) Ax- 
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Applying this to the tq, Ax- sa case shows that £^ = 1/^/1 — ^ so that finally 

2—2 I A— |2 7p2 ( 1 '^'^ \ 2—2 2—2 

cr„-|Ax„| |l--^)cr =ct 

Hence equation (|3.12() applies to the drift /time pair Ax^, r„ and AaJp « 0, tq hence 
Ax,\ ^ Mx7o 

1 /I — w \ / Ax-Q 



c' 

Moreover equation (|3.12|l relates the drift/time pair Ax, t and Axo ~ 0, tq so 
again 

^Ax\ _ /Ax-„^ 
r ; - 1,£^K) j ^ TO 

_ 1 / 1 -v'\ fAx^: 



Finally from these two equations it is clear that 

AxA _ (A{v) B{v)\ (A{v') B{v')y^ fAx^ 
r„ ; - [f{v) E{v)) \F{v') E{v')) \ r 

which yields 1)3. so the argument is complete. □ 

Remark 3.4. The analogy with relativity extends further than just the Lorentz 
transformation above. For example, for the kinetic energy of the main particle it is 
possible to write E \Ht] ~ \Mtv'^ « ^Mv"^ ignoring the osmotic term, the eflFect 
of a finite mass ratio and the diffusion term in (j2.15|l . Expressing the differential 
|Axt| in terms of To 

Em^M^-^^MU--!^^ 



and recalling the definition of the relativistic energy i?,. = M(? j ^\ — jc? this 
equation can be rewritten as 



M \Ax-f Mc^ 



which yields 



= mJ-^^ + Mc^ = + Mc^ (3.15) 



with Mr = , ^ „ the relativistic mass. 
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The remark allows an interesting interpretation of the various terms E [Hk], c 
and the diffusion term e/(a^A/r). Clearly the relativistic mass Mr in this context 
is the mass measured as a function of correlation while the kinetic energy E [Hk] 
refers to the total kinetic energy of the associated main particle and is in Relativity 
Theory referred to as the relativistic momentum. The second term on the righthand 
side of (|3.15(l is the rest energy of the main particle and is in the present context 
the result of the diffusion energy. In fact, using H3.6|l it is clear that 



interpreting the relativistic mass in terms of the mass ratio 7, the average inter- 
collision time To for the main particle at rest in the heatbath and the variance per 
mass ratio e (with units of action). 

The most important conclusion from these observation is that it is possible 
to find a simple relationship between the inter-particle collision time and for- 
ward/backward (main) particle velocity correlation in the form of the Lorentz 
transformation, see Lanczos [17] for some more detail. The relativity analogy 
can be pushed further and the relativistic momentum and the rest mass energy 
of the main particle can be interpreted in terms of the time to collision, mass 
ratio and variance per unit mass e. Unfortunately, the compounded correlation 
pf, = (1 — 2a^) Pa+zA-z destroys the martingale property for the coordinate pro- 
cess of the main particle and this in turn affects the correlation structure of the 
heatbath particle. Theorem (|2.8|) shows the heatbath backward and forward veloc- 
ity are correlated already so the superimposed correlation will have an additional 
effect as shown below. 

The remaining part of this Section will briefly discourse on a correlation model 
for the impulses ■^aA'^z and -^aA'z. A convenient route is to assume that 



(3.16) 



where A^-z^ , ArZ~' are independent Gaussian increments, tTa = aa{x,t), Uo — 
ao{x,t) are n x n matrices and Aq G M" and Aq G R" independent processes 
with E [AaA^] = E [AoA^] = 2//r. Proper conditions on the drift and vari- 
ance terms will not be specified here. This representation is motivated by the fact 
that by this construction A+a;(t, /3) — A~x(t,/3) does not depend on aaAa and 
A+a;(t,/?) -I- A~x{t,P) becomes independent of CToAq. 

The form of the correlation can now be summarized in a straightforward calcu- 
lation as follows. 

Proposition 3.5. If the correlation structure of the main particle is represented by 
equation H3.16|l then the random variables Aa and Aq are part of both the past and 
the future steps of the process or rather relate the past to the future. The correlation 
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structure of the heatbath particle is then represented by 

"^-^^ ^ ' cos{0) 1 

1 - cos(6') 





(3.17) 



7sin(0) 

where A,.z+, ArZ~ ,aaix,t),ao{x,t), Aq G M" and G M" are defined above. 

Proof. The first remark in the proposition is seen from the fact that equation (|3.16() 
imphes that 

„ ^ f A+z A"z\ /A+z A::z~ 

2craAQ = cr 



. " X . ! ( (3.18) 

/A+z A-z\ /A+z A-z^ ^ ^ 



The matrix in the proposition is the resuh of inverting equation l)2.9b|l so (|3.17(l 
follows from the fact that the Ao G K" and Aa G M" terms are eigenvectors. In 
fact, 

1 /-C0S(6I) 1 \ faaAa\ _ (cTaAa 



7Sin(6') V 1 -C0s(6l)y \aaAa) \(TaAa 

1 cos(0) 1 \ / aoA, \ _ I f cToAo 



jsm{0) \ 1 -cos(6l)y \-aoAoJ 72 \-aoA, 

since = 1 and "'"^^^^g^^'' = 1/7^- This completes the Proposition. □ 

To show the effect of the ctqAo and <TaAa terms take for example n — 1, then 
equating the variances on both sides of H3.16|l 

—a'^var (A^z) = cr'^var (A~z) = 

2.2 (3.19) 

alvar {Aa) + alvar (A^) + 



T 



Since A^ and Aq are independent a'^E ^^^r^ = alvar {Aa) — a^var {Ao) so 

that the correlation between between A+z and A+z reduces to 

alvar {Ag) crgmr (A^) 

PA+2A-2 = (,^.z!Uj 

cr^uar (A^) + alvar {Ao) + ^ 

If it is assumed that var (A^) = 2/t = war (Aq) = 2/t then (|3.2Q|I reduces further 
to 

a^ ^2 

The reason for this construction now becomes clear. If ctq « and (Tq >> f r then 
Pa+zA-z — > — 1 while Pa+zA-z ^ 1 if ctq w and da >> cTr- The correlation in the 
driving factors in H3.16|l relates the future to the past through the collision time t. 
This introduces a form of auto-correlation for the random difference process. For a 
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positive or negative correlation the process x {t, f3) ,t > can not be a martingale 
and a two dimensional process must be introduced. 

This Section relates correlation in the motion of the main particle to the inter- 
particle collision time via the Lorentz equation to satisfy a Minkowski invariant 
and suggest that the rest mass energy is the result of the energy embedded in 
the diffusion energy. Proposition H3.5(l shows how to decompose the backward and 
forward velocities into a perfect correlation part and a martingale process. These 
results were predicated on the interaction prescription (|2.()l) which is exact in the 
one-dimensional case. In higher dimensions the equation applies to the center of 
mass line projection as remark (|2.4() describes. The energy conservation arguments 
in the last two Sections still apply however another level of complexity will be 
required to describe the non-simple collision. This will be addressed in the following 
Section. 



This Section returns to the case where the matrix fl in equation H2.6|l char- 
acterizes all elastic interactions incorporating the random anti-symmetric matrix 
Z = Z{U). In this case the total energy defined in Theorem H2.5(l depends on 
the statistical characteristics of Z and Proposition l|2.1()|l is derived in the presence 
of this collision scattering matrix. Finally the conservation condition for the full 
collision is derived and an example is presented that combines this result with an 
electromagnetic field type Hamiltonian. The Section concludes with showing that 
it is possible to subsume the collision scattering matrix Z into the heatbath mak- 
ing all the results from Section 2 and 3 applicable for an altered heatbath with a 
different statistical structure. 

The first step is to obtain the total kinetic energy expression in Theorem H2.5|l as 
a function of the pre- and post collision velocities V2 , vi of the main particle using 
equations H2.2|l and H2.5|) . The form of the total kinetic energy is introduced in the 
following result. 

Theorem 4.1. As in Theorem H2.5(l . let the momentum of the main particle and 
interacting particle be presented as pi = Mvi (post- collision p2 = MV2) where 
Pi,P2,Vi,V2 G M" and qi — Mwi (post- collision (72 — MW2) with qi,q2,wi,W2 S 

1 2 2 1 2 2 

R". Then the total kinetic energy Tik — i^iM \v2\ -\-m\w2\ ) = ^{M \vi\ -t-TO|ti;i| ) 
is related to the pre - and post collision momenta of the main particle as follows 



4. Non-Simple Collisions / Scattering 



8Hfe 



Mt 




(4.1) 



(A+t; - ZA-y)'^ {A+v - ZA'v) 
+ ^A'v'^A-v + ^A-v'^ZZ'^A-y, 



where 



A+U = V2+V1, 
A^V = V2 — Vi. 



(4.2) 
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Here Z = Z{U) is the anti- symmetric matrix (Z'^ + Z = 0) such that Z = I — 2{I + 
U)^^ = / — 7sin(0)Q^^ where the unitary matrix U and Q are defined in equation 
(|2.5|) above . 

Proof. The proof can be found in Appendix B as well. □ 

As mentioned above the random collision scattering matrix Z = Z {U) represents 
the collection of all possible elastic collision transitions and therefore contains the 
center of mass line information and the angle of impact between vi and wi . The 
matrix will be different from one collision to the next and has a statistical mean 
E[Z] = Z and variance E [ZZ"^] which may be a function of the collision position 
x{t). To deduce the expectation of the energy term Hk use again H1.3|l with the 
definition 1)4. l|l to derive 

8E [Hk] ^A+v'^A+v] - 2E [A+v'^ZA-v+] 



Mt 



„2 



, E[A-v'^A-v] + (^^]e[A-v^ZZ^A-v] , ^ 

Y ^ ^ V 7W (4.3) 

=E [A+v'^A+v] - 2E [A+v'^ZA-y] 
+ ^E [A'v^T'A-v] , 

where E[Z] = Z and = / + (1 + 'j'^)E {ZZ'^). This is the result of taking the 
expectations over the random matrix Z first and then rearranging the expression. 
Alternatively, this expression can be written as 

SE [Hk] 



Mt 



--E 



{A+v - ZA-vf {A+v - ZA- 



+ ^E[A-v^T^-^^ZZ)A- 



T 

Notice that the matrix in the last terms is positive definite since — j'^ZZ = I + 
ZZT + ^^ (Yzt - -f^'ZZT^ = I + zW + j^var {ZZ^) with the obvious definition 

for var {ZZ'^) = E [ZZ^] ~ E[Z]E [Z^] > 0. 

Again it is assumed that the change in expected energy equals the change of an 
appropriate potential $p so that j^^E \Ht] = jj^E [Hk + = 0. The following 
Proposition shows the form of the total energy H4.3|l as a function of the backward 
and forward velocity. 

Proposition 4.2. Reducing the expectations in H4.3fl the total energy in can he 
expressed as 

T 




( 



b++b- 
2 
T 



z 



—E[Tr{T^)] + —E[% 



(4.4) 

1 

2t ' 277^" ■ " • Mt 
Proof. This expression can be easily derived from substituting l|4.2(l into 1)4. 3|l . Let 
zi = ^ ^ ^ ^^'^^ ^""^^ ^'^'^ independent and 
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normally distributed. Hence 



4n 



E [zf zi] = E [zjz2] = — / 



E 



A- 



T2 



T 

A+z 

T2 



A- 



An 

T 



E [zlz2] = 0. 



Moreover, it is easy to see that E [z'[Z'^Z2] — and E [zjz2 + (l + 7^) zJZ'^Zz2] — 
^E\Tr{V-)\. □ 

The pre - and post -collision velocities for the main and heatbath particles are 
linearly related via equation 12.2|l for the matrices P^Q^V and G defined in H2.5|l . 
Now the equivalent of Theorem (|2.12|) is introduced to show the conditions for 
maintaining a constant total energy. 

Theorem 4.3. Let p — e "'-^ — e ^ and introduce the sufficiently smooth func- 
tions A — A{x,t)^S = S{x,t),x £ R" and constants (5,^ to express the backward 
and forward drifts 6+ b'^{x, t),b^ = 6^ (x, t),x £ M", t > Q as follows 

6+ = e (yS -A)+ jSVR, 
5- = ^ {WS -A)- -fSWR. 

Assume that the potential $p satisfies the following property 



j^E[%]^e[{VS~A).{V^ + ^A 



E 



A. 



(4.5) 



with A — ^ using Einstein's notation of summing all like indices. LetTCk be defined 
as in equation (14.4(1 with Z a random matrix such that = / + (1 + 7^)i? [ZZ'^'\ 
and E [Z] = Z . Then the total energy "Ht = Ti-k + for the potential in $p in 
(|4.4(l is conserved if 

d_ f E[nk] + E[%] 

dt 



1 /^C^ IVS* - A |2 - 2e(57(VS' - A^Z^BX d 1 



-Jt2^ 



(5, 



Ar. 



dx 



(4.6) 



J p{S.x^~Ap) (e-^) Apdx + EiZ) 
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whe 



+ I / P (5.. - \) ((^'^7 {S., - A,) dx (4.7) 

- ^(57 / p iJySx^t - ZjkRxt + {Sxj - Aj) ZjkRxk^ dx, 



and where Tr (F^) denotes the Trace of the matrix . Here 

Proof. For a proof consult Appendix C. □ 

This representation exhibits three sizeable problems finding solutions for the time 
invariance of (|4.4(l . First of all there is the fact that potentials $p — ^p{x, t) do not 
typically admit property like l|4.5|l as the time derivative of the potential introduces 

Obviously time independent potentials satisfy property (|4.5(l 



terms like E -g^ 

and as can be seen below Maxwellian type fields have this property as well. In fact 
if the Z term can be ignored for a moment then the following can be shown. 

Theorem 4.4. Assume that Z = and let ip — i}){x,t) — e ^ with 
p(x,t) — \%l){x,t)\ where 

6+ = ^ (VS' -A)+ -fSVR, 
6" = ^ {WS - A) - -fdWR, 

and let E and B satisfy the (magnetically sourceless) Maxwell equations 
S/.E^p{x,t), 



S/xB-^^^E = i(^^^^]p{x,t), 
V.B = 0, 



(4.8) 



so that E = -V0 - £,-§^A and B = V x A. Assume that $p = \Ef + \Bf then 



'dt 

d fE[nk]+E[%] 



dt \ Mt 

if and only if 



= 0, 



ixA^-:^{x^ -^Af^ + c^{x,t)i;, (4.9) 



with X = MtT] = Mto-'^M = + 7) ^ '^'^^ 5 = ^ = l/Mj 
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Proof. The point of the proof is that for equations (|4.8(l it is true that 



E 



(v0 + ^i 



In other words the time change of the energy of the field in (|4.8f) satisfies (|4.5ll and 
this is combined with proposition (jC.ip and equation ljC.5|) from Appendix C. □ 

The second issue is the fact that the terms ^^RxjTjf.Rx^.,kj,k — and 



,fc = 1,, 



, in Theorem (|4.5|l above depend on the random Z matrix 



T 

which acts here as an arbitrary scahng factor. If the matrix is diagonal is is 
possible to scale the solution derived for the case where Z is state independent. In 
fact the example below shows that a simple scaling applied to the factors in the 
wave function allows for a solution. 

4.1. The two-step Scattering Matrix. This example derives a solution to equa- 
tion (|4.t)|l by rescaling the wave function as follows. Assume that n = 2, ^ = 
A = (f> = and define the anti-symmetric matrices Z{i'), Z{~v) as follows 



Z{y) 



V 

-V 



.Z{-v) 



-V 
V 



Now let the probabilities -piy) and p{—i^) be such that 
E[Z] = Z = p{y)Z(y) +p{-y)Z{-y), 
E[ZZ'^] =p{v)Z{u)Z'^{v) +p{-v)Z{~v)Z^{-v) = i^^I, 

so that = 0, Tr (j^^^ = with Z state independent (not a function of x{t)). 

Then = / (l + (l + 7^)i^^) = la^ and the solution that preserves the energy 
equation H4.6|l reduces to 



d 
di 



E[n] + E 

Mt 



Sx„ 



^St + ^ |V5| 



|Vi?|- 



dx 



(4.10) 



because the second term in equation (|4.fi|l vanishes 
If / pSx, ((C^75.,^,.).,.J dx 



dx = 0, 
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and so do the last two terms. Hence the solution to H4.9(l is given by Proposition 
(|C.3|I below as 



2Mt 



(4.11) 

where — iplx^t) = e ^ with Xi' — ^Itv/'^i^j C — ^/Mt and (5 = 

1/ [a^MT)- Proposition H4.10|l in Appendix C has some more details on the deriva- 
tion. 

This highlights the second feature of equation (|4.6|l which is that the effective 
osmotic term typically scales up due to the variance of the random matrix Z . This 
scaling never disappears if E [ZZ'^'^ > however equation (|4.6|) can still be reduced 
to the case of the simple collision in Theorem (|2.12() as example (|4.1(l shows. There 
is no real quantum mechanical analogy to this except that some of the terms in 
equation (j4.10|l are similar to the terms in the Bopp-Haag Hamiltonian where the 
additional derivatives are introduced to incorporate spin states, see for instance 
Nelson fH. 

The third question on equation H4.6|l is the effect of the Z terms. The example 
above shows that these terms disappear for the case where Z does not depend on 
the coordinate system. This is entirely due to the fact that Z is anti-symmetric and 
their contribution to equation (|4.6|l remain conveniently zero even if the matrix Z 
depends on time t. As the mean collision scattering matrix relates to the average 
center of mass line and average angle of collision this quantity is unlikely to be 
dependent on the collision coordinate x(t) so this assumption is not unreasonable. 

Obviously the presence of the Z matrix changes both the correlation between V2 
and vi and simultaneously affects the correlation structure of the heatbath. If for 
instance the main particle path is a martingale for Z = then "turning on" the 
Z will create a correlation. If the collision scattering matrix Z > and the main 
particle path is a martingale then the correlation structure of the heatbath must 
change from the case that Z = 0. The following Proposition generalizes Theorem 
(|2.8|l and calculates the correlation of the heatbath particles for the latter case. 

Proposition 4.5. Assume that the backward and forward velocities of the main 
particle are uncorrelated. Then the correlation matrix for the colliding heatbath 
particle given a realization of the random matrix Z looks like 



E 











f) 






\ T2 






2 / 


1 



^1 - 2a-) 



- (1 - 2a-] 
1 



whe 



2a- 



r7 sin (9) 



E [ZZ^] 



E IZZ^] 



id where 



17 = (1 - cos(6l)/) Z + E [ZZ^] . 
Proof. For a straightforward calculation see Appendix D. 



□ 



ENERGY AND MOMENTUM CONSERVATION FOR STOCHASTIC PROCESSES 43 



Remark 4.6. Though Proposition (|4.5(l specifies how the heatbath particle field 
must behave to guarantee that the main particle moves in a Markovian fashion the 
distribution of x{t) has now become complicated. Clearly the distribution for V2, W2 
must be convolved with the distribution for vi, wi and Z = Z{U) and is therefore 
not readily calculated. 

The final question in this Section is whether there is a " momentum term" version 
of Theorem and whether the complexity of equation (|4.6(l can be reduced. 

Interestingly this can indeed be achieved in a straightforward manner but some 
changes in assumptions will be required. The approach is to absorb the scattering 
matrix Z into the heatbath and then show that the results from Section 2 and 3 
apply for the transformed heatbath. 

To implement this approach the A+i; and ISTv must be rewritten in a more 
convenient form. From the proof of Proposition (|4.5(l in Appendix ^ it becomes 
clear that 

W2\ _ 1 /- C0S(6I)/ -Z I + Z \ (V2 



^wij 7sin(6i) V - cos{9)I + Z J \vi 

This can be simplified by writing 

W2\ _ (Vvi + Gwi\ _ fV G\ fvi 



and then use equation ljB.8|) and (jB.9|l from the Appendix to obtain 



1 {V G\ (1 -I \ fA+v 



W2 

ywij^2{0 I J {I {2Q-^-l)J {A-v 

hence with Vr = (2(5~^ — /) this reduces to 

W2\ _}_(V G\ (I -I\ (A+v 
Wi)^2\0 l)\I Vr)\A-v 

I (V + G ~V + GVt\ (A+v 



2\ I Vr J [a 

I (I Vr~ (A+v 
2\I Vt J {a-v 



Inverting this relation yields 

= 2 



(A+v 
\A-v 




which means 



A+v = 2wi + 7^yrA"w, 



A v — — 7^A w. 

Now -f'^Vv = / - (l + 7^) Z so then 

A+v = A+w - (l + 7^) ZA~w, 
A~v = — 7^A~w, 



(4.12) 
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which can also be written as 



A+v = A+w + 2>V, 
A^v — —j'^A^w. 



(4.13a) 
(4.13b) 



Here mW is the additional momentum such that — W = Z{w2 — u'l) 



1+7 



ZA' 



From this the following two Theorems are easily shown. 

Theorem 4.7. The additional momentum defined in equation H4.13a|l above is 
perpendicular to the A^w, in other words 



^ - ^ + ^—A-uFz^A-w = 0, 



A-yj = 



and so if E 



E 



\W2\ 



then 



E 



= E 



\W2+WY 



(4.14a) 



(4.14b) 



Let wt.i = wi + yy, wt.2 = wi + W and define A+v = wt.2 + wt.i, A^v 
wt.2 — u)T,i = W2 — Wi then equation H4.13a|) reduces to 



A+v = A+w, 
A~v — — 7^AyW 



henc 



where = E 



\wt,2Y 



1 

2^ 



E 





A+v 


to 


1 


A^v 


2' 




2 






2 





4, 



(4.14c) 
(4.14d) 



(4.14e) 



Proof. To prove assertion 14.14a|) consider (|4.13a|) and multiply with (j4.13b|) so that 

A~v'^A+v = A~v^ {A+w + 2W) 

= --f'^A~w'^A+w - 7^A~w;^>V, 



or 



\v2\" - \vi\" + \w2\" - hiT = -7^A-w^W 



However the left hand side equals |i'2|'^ + 7^ 



f-{\vi\' + l^\wi\') =-Hk/M- 

Tik/M = by energy conservation hence 7^A~w-^yV' = and H4.14a|) is proved. 
Assertion (j4.14b|l now follows since 



E 



E 



\W2+WY 



= E 



- E 



\W2\ 



2E 



{wi — W2) W 



-2E [A-nFw] = 0, 
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and equation (|4.14e|l follows directly from (|4.14c|l and (|4.14d|l since 

= A^+w^A^+w + A^^w^A'^-w 



1 T|2 









This concludes the proof. □ 

In other words, the momentum constraint on the forward and backward drift 6+ 
and 5~ with Z ^ with the heatbath particles moving at average speed c can be 
transposed into the case where Z = with the main particle in a heatbath where 
the particles move at average speed ct- Here again — mW = ^^^mZ{'W2 — wi) = 
— '^'^2 mZA^'w is the additional momentum. Moreover, due to l)4.14b|l and H4.14e|l 
the momentum conservation requirement for wt,! and it;T.2 as described in Theorem 
()3.3|) is equivalent to a requirement on wi and 'W2- The sole difference between the 
original and the transposed case is that the correlation between wt,i and wt,2 
is different from the correlation between wi and W2 due to the fact that W is a 
(random) function of Wi and W2- The exact expression is calculated in Proposition 
(14.51) for the case where the main particle follows a Markovian path. Any changes in 
correlation going from wi, W2 to wt,!, wt.2 suggest that the average inter-collision 
time is affected as detailed in Section 3. 

The next result completes the identification between a heatbath W2 , wi in which 
Z and a heatbath wt,2, wt,i in which Z = 0. While (|4.14e(l is the transposed 



heatbath equivalent of the momentum constraint l)2.7b|l it is not immediately ob- 
vious that H2.7a|) has an equivalent as well. The next Theorem shows that this is 
the case and as an addendum calculates the original energy in terms of the original 
heatbath terms. 

Theorem 4.8. Let the main particle diffuse in a heatbath where wt,i = wi + W , 

WT,2 ^ W2 + W with -W = ^^Z(w2 - wi) = -i^ZA-w. Then the total 
kinetic energy TLt ,k equals 



= A+v'^A+v + ^A-v'^A-v 
M 72 



\ \V2\^ + \ |WT,2|^ 



(4.15a) 



- A+v^A+v + -,A-v'^A-v 
M 74 (4.15b) 

+ 7^ A+w^A+w -I- -f'^A~'uFA~w 

Proof. No calculations are required to show assertion 14.15a|) as the definition of 
yy in H4.12|l and (|4.13a() shows that 

^V2\ _ 7 (C0s{e) 1 \(WT,2\ ^^ -^g^ 



vi ) sin(6i) V 1 cos(6') ) \wt,i J ' 
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using (|D.3(I and (|D.4() . From this it is easy to work back and derive 



^cos(6') 7sin(6') \ f vi 
cos{9) j \wt,-i 



V2 
WT,2 

which imphes in turn that (|2.7a|l holds with wt,i, u't.2 replacing wi and W2- Hence 

2 





1 


V2 + Vi 


2 


1 


V2 - Vi 


Mt 


~ 2 


2 


H 


' 2^ 


2 



with Ht k = + T I^T i|^ = ^ 1^2 1^ + ^ |wt 2!^ but this is exactly assertion 

Equation (j4.15b|) can be shown from straightforward calculation. Equation 
(IB. 1311 shows that 



Mt 



A+w^A+w - 7''A"w;^ A^w 



2A+?i;^ZA"w + (1 + 7^) A^w^Z^ZA 



w, 



but from (|4.12|) it follows that 



Finally then 



A+v'^A+v + —A^v^A~v 



A+w^A+w + 7^A~w^'A"'u; 



+ (1 + 7^) 



~2A+uF ZA-w 
+ (1 + 72) A+w'^ZZ^A+w 



A+w^A+w + -^'^A^vF A~w 



+ (1+7^) '\,2^-...T 



- A+w^A+w 
j'^A~w'^A'~w 



A+w^A+w + ^A~t>^A"i; 



- j^A+w'^A+i 



M 

— j'^ A^ w'^ A^ w , 



which concludes the proof. 



□ 



5. Conclusions 

The purpose of this paper is to revisit the classical diffusion problem of a main 
particle moving through a heatbath propelled by elastic collisions and introduce in- 
teraction energy and momentum considerations. The main particle path is modeled 
as moving linearly from one random collision to the next with random inter-collision 
times Tj,j > 0, represented as second order Gamma distributions. This pre-limit 
microscopic construction models the motion of the main particle by collision posi- 
tions x{tj),j > 0, at stopping times tj, j > 0, with tj+i — tj ~ tj, j — 0, 1, .... The 
derived Hnear interpolated path x {t, (3) ,t > constitutes a "best estimate" of the 
main particle position. An important model assumption here is the choice of the 
second (or higher) order gamma distributions as the particle inter-collision times 
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Tj,j > 0. If the distance traveled between the colhsions x{tj) and x(tj+i) is mod- 
eled as b'^{x,t)Tj and a (subordinate) Gaussian contribution then the inter-collision 
momentum and energy of the main particle is finite with probability one. 

Section 1 shows that if the mean inter-collision times decreases the main par- 
ticle path approaches the strong solution to the continuous stochastic differential 
equation with drift b^{x,t). The drift and the variance term must satisfy standard 
requirements which guarantee the existence of a strong unique solution but one ad- 
ditional time growth condition on the drift was introduced to control the increasing 
multitude of "collision" path contributions. Some constraint must be applied to 
the time-variability of the drift because otherwise there is no guarantee that the 
drift specified at the collision points does not deviate too much from the continuous 
drift function. 

Section 2 follows the consequences of the fact that the collisions with the heatbath 
particle are elastic and introduces the canonical solution to the collision energy and 
momentum conservation constraint. The pre- and post collision velocities of the 
colliding main and heatbath are linearly related via a matrix containing a random 
collision scattering matrix Z = Z[U) which carries the center of mass line and 
impact angle collision information. For simple collisions for which Z = 0,U = I 
the canonical solution implies two constraints on the motion of the main particle. 
The first one relating the total kinetic energy (main and impacting heatbath particle 
combined) to the pre- and post-collision motion of the main particle. The second 
relationship has been referred to as the "momentum constraint" and looks more 
like a velocity requirement. If the total kinetic energy along the path of the main 
particle is not constant (on average) then some energy is being transferred between 
the main particle and the heatbath. 

The main result in this Section shows that if there is no energy leakage (on av- 
erage) between the main particle and the heatbath hence if the total kinetic energy 
plus possible potential is conserved then the probability distribution of the position 
of the main particle must be derived from Schrodinger's equation. Planck's con- 
stant is then replaced by a variance per unit of mass term e and further depends 
on 7 the mass ratio. The total energy functional may contain only certain suit- 
able potentials satisfying an energy conservation property. The derivation relies 
heavily on the stochastic mechanics results and on convergence of the collision path 
representation to a suitable stochastic process. 

The important aspect of this derivation is that a combined energy constraint for a 
particle diffusing through a heatbath is a purely classical problem and Schrodinger's 
equation is invoked to prevent energy exchange between the main particle and the 
heatbath. The analogy with quantum mechanics however is not perfect as for in- 
stance the forward and backward drifts explode when the mass ratio 7 becomes 
small. The present derivation moreover only allows certain appropriate potentials 
though this includes all time-independent potentials and the electromagnetic po- 
tentials see Section 2 and 4. Apart from the presence of the mass ratio 7 there is 
an important conceptual difference. The total energy functional that is conserved 
by the presence of the Schrodingcr wave function is the total kinetic energy of both 
the main and heatbath particle. Quantum mechanics associates the wave function 
only with the "main" particle and the energy eigenvalues of the wave function sup- 
posedly arc the energy states of the main particle alone. In this paper the energy 
states obtained from the wave function relate to the combined main particle and 
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heatbath particle energy. Once the wave function solution has been obtained the 
main particle and heatbath particle energy contributions still need to be separated. 

Yet another difference between the present formulation and the Theory of Quan- 
tum Mechanics is the curious aspect noted in the Gaussian Wave Packet example 
(|2.2|) where the statistical characteristics of both the main and the heatbath parti- 
cles are affected by the energy constraint. The kinetic and energy dispersion terms 
that are found in the energy of the main particle can also be found in the heatbath 
energies so the heatbath is affected by the presence of the main particle in a time 
dependent manner. This effect decreases if the mass ratio 7 becomes smaller but 
the model can not be structurally amended. The results in Section 3 suggest in fact 
that a different approach may be required which incorporates correlation between 
the main particle pre- and post collision velocities. 

Intuition suggests that the diffusion/quantal effect becomes noticeable if the dif- 
fusion per mass term is significant in comparison to the size of the drift term or the 
energy dispersion. The example of the Gaussian Wave particle in Section 2 shows 
that the energy dispersion and the average main particle momentum increase the 
combined kinetic energy. This suggests that a high temperature heatbath environ- 
ment that is not not overly dense is dominated by its energy dispersion and mean 
particle drift. On the other hand if the main particle moves through a relatively 
narrow energy band with a relatively small kinetic motion then the diffusion term 
constitutes the larger part of the total kinetic energy. Ultimately if the heatbath 
is very dense then the diffusion term will become the dominant energy provider. 
A very dense heatbath environment eventually forces the dominant motion of the 
main particle to be entirely diffusive. 

The example also shows that the heatbath looks like a reflection of the main 
particle. The main particle carries almost all the momentum and dispersion energy 
with little diffusion energy while the heatbath has a very large diffusion term and 
very little kinetic energy. If the pre- and post collision velocities of the main particle 
are uncorrelated as one would expect for a Markovian path then the pre- and post 
collision velocities for the heatbath particles will be correlated through the collision 
point. Both the main and heatbath particle kinetic energies are time dependent but 
the first one increases in size while the heatbath energy decreases proportionally. 
For a very small mass ratio 7 the time dependence almost disappears and the 
heatbath starts to behave as if it has only one velocity which is reflected by the 
collision. 

The momentum constraint established in Section 2 shows that a similar sym- 
metric quadratic expression employing the forward and backward velocities of the 
main particle can be directly related to the average velocity of the incident parti- 
cle. The average is calculated as the arithmetic average of forward and backward 
heatbath particle velocity. The contribution of Section 3 is to show that if the 
heatbath particles are in energetic equilibrium with the main particle then the drift 
of the main particle and the correlation between the forward and backward veloc- 
ities of the colliding heatbath particles must depend on the average inter-collision 
time. The condition is identical to the geometrical Minskowski invariant employed 
in Special Relativity and it is shown that the invariant can be satisfied by applying 
the Lorentz transformation to the average collision time and the squared distance 
traveled. 
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The analogy with Relativity Theory can be pushed further to suggest that the 
energy rest mass of the main particle equals its diffusive energy. In fact the rest 
mass can be expressed in terms of the mass ratio 7, the average inter-collision 
time To for the main particle at rest in the heatbath and the variance per mass 
ratio e. This argument only employs the Lorentz transformation as a means of 
generating a solution to the Minkowski invariant and is not necessarily the only 
solution which balances the mean inter-collision time and the squared distance 
traveled. By comparison in Relativity Theory the homogeneity of space and the 
constancy of the speed of light in all directions leads to a unique solution. 

The third set of results in section 4 focusses on the " non-simple" solution to the 
elastic collisions to include the random collision scattering matrix Z . The path of 
the main particle now becomes a " conditionally" Gaussian process in the sense that 
the main particle path process remains a Gaussian process given the realizations of 
the random collision scattering matrix Z . In general however the additional random 
matrix Z alters the main particle coordinate distribution. The stochastic dynamics 
for the random matrix involve the center of mass line distribution and depend on 
the dimensions and physical setting rather than on the motion of the main particle. 
It is therefore reasonable to assume that the collision scattering matrix Z does not 
depend on the pre- and post velocities of the main particle. 

One of the most obvious effects of the random collision scattering matrix is the 
arbitrary change in correlation between the pre- and post collision heatbath veloc- 
ities causing in turn a correlation between the pre- and post-collision velocity of 
the main particle. This destroys the Markovian property of the main particle path 
process so that additional variables must be introduced to model its motion. None 
of the results in Section 2 are then applicable because the fundamental Markovian 
result relating the backward and forward drift difference and main particle position 
probability density is not valid. However Section 3 demonstrates how the pre- and 
post-collision velocities correlation of the main particle process relates to the mean 
inter-particle collision time which can then be captured by the Lorentz transforma- 
tion. Unfortunately there are no results that describe how the probability density 
of the main particle position can be calculated in the presence of correlation. 

Theorem (|4.3ll in Section 4 presents the main result showing that the conservation 
of total energy depends on the mean scattering matrix Z and expected covariance 
matrix E [ZZ'^'\ . A full solution for the probability density has not been derived 
but for the case where Z = it was shown again that the probability density 
must be obtained employing Schrodinger's equation as long as the potential is 
time independent or satisfies a Maxwell type set of equations. Another example 
is presented for the case where Z is a function of time only while E [ZZ'^] is a 
constant diagonal matrix. Then it is possible to obtain the probability for the main 
particle position in the form of a wave function for the Z = case but with some 
of the weighting parameters altered. In comparison to the Z = wave function 
this solution looks as if the mass weighting and the diffusion per unit of mass have 
changed to account for the collision scattering matrix. 

A very useful result from this Section is that it is possible to absorb the collision 
scattering matrix Z into the heatbath. Specifically it is possible to represent the 
case with a non-zero scattering matrix Z ^ in a wi, W2 heatbath with the case 
where Z = with a w-f.i — wi + W, wt,2 — W2 + W heatbath. The additional 
momentum mW is proportional to the collision scattering matrix Z and acts as 
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a straight increase in the energy in the heatbath while it is orthogonal to the 
momentum exchange. However the correlation structure for wt,i, wt,2 is different 
from the correlation matrix for wi and 102- Therefore the path for a diffusing main 
particle can not be Markovian in both the wi, W2 heatbath and the wt,i, wt,2 
heatbath simultaneously see Proposition H4.5|l . 

The results in Section 2 do not apply when the main particle path is not a 
Markovian process however the results in Section 3 take correlation into account 
as long as the correlation is homogeneous. In fact Theorem Ij^.tfl establishes a 
relationship between the mean particle speed = \wi\^ = |it;2|^, the compound 
correlation py and the mean inter-particle collision time r. As a result of (|4.14b|) 
this exact relationship must hold for for some Cy = |u't.i| |'i«T,2| , p'v and t' as 
well. The results of Section 2 may still provide a good Markovian approximation if 
the Wt,i, wt,2 heatbath renders the pre- and post-collision velocities of the main 
particle independent. 

From section 3 it is clear that the motion for a main particle in the presence 
of correlation between pre- and post-collision velocities must satisfy a relativistic 
invariance hence extrapolating from results in Section 2 and Section 4 a proper 
distribution for the main particle position is likely to satisfy a type of Klein-Gordon 
equation. This line of research should pursue the ideas of Serva or Guerra ^1] 
and will be investigated in Part II of this paper. 
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Appendix A. 

Proof of Theorem (j2.5(l Equations H2.7a(l can of course be verified by direct 
substitution however the following simple argument is more intuitive. Using pi = 
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Mvi, p2 = Mv2, qi = Mwi and q2 = Mw2 it follows from equation 1)2.6(1 that 

92 + 91 = Pi7 sin(0) - qi (cos(0) - 1) 
= Pi7 sin(6') +qi2 sin (9/2^ 
= 2pi7sin (9/2) cos (9/2) + 2qi sin {9/2f , 



92 + 91 



= sin (9/2) (pi7 cos (61/2) + qi sin (61/2)) . 



Similarly using momentum conservation P2 ^ Pi — ^ (92 ^ 9i ) ^nd equation (|2.t)|) 
results in 

7 — ^ — = -pi(l-cos(6l))-qi-sin(6') 

= Pi7 (sin (61/2))^ - gi sin (6'/2) cos (6'/2) 
= sin (61/2) (7pi sin (61/2) - qi cos (61/2)) , 

and adding the squares results in 



92+91 



92 - 9i 



7^ sin (0/2)2 bilV 



I9i| 



Substituting 



into this yields 



hence 



2Mn 



2Mj* 



k-f'sin{9/2y = —^Hk 
1 + 7 



92 - 91 = Pi -P2, 

92 + 91 = 7^(Pi +P2), 



P2 +Pl 



P2 -Pl 



P2 +Pl 


' 1 
+ — 


P2 -Pl 


2 


7 


2 



21/7" 

1+7" 
2M 



1 + 7- 



which covers the first two expressions in Theorem (|2.5ll . Now dividing this equation 
by yields 

2 2 

M(l + 72)^'- = Mj:^'^' 

proving equation H2.7a|l . 

Direct substitution can be employed again to obtain equation l|2.7bp but this 
equation can be derived more organically from the pre- and post-collision total 
kinetic energy and equation H2.7a|l . 



V2 + Vi 


' 1 


V2 - Vi 


2 


+ — 

T 


2 
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So for a direct derivation write 



^fc = ^ 11^21^ + 7^ |W'2| 



Af(l + 72) 



W2 + I'l 


^ 1 


W2 - Vi 


■) 


2 


272 


2 





(A.la) 
(A.lb) 
(A.lc) 



so that 



{ 


V2 + Vi 


' 1 


V2 - Vi 




2 


^ 272 


2 



^(h|VK|') + M72K|V^'^^' 



M , 



\W2\^ - \wi\' 



2 " ' ' ' ' ' ' ' 2 

which is achieved by adding IjA.lap and IjA.lcp and equating that to twice the 
middle term (|A.lb|) . 

Carrying the V2, vi terms to the left hand side then yields 



M 9 9 

Y(k2|VKr) 



M(l + 72) 



(l^2| 



1 



2u9 vi 1 



M 9 9 

Y(k2|VKr) 

2 I |2 , l\ ,2 I |2 

' ' ' 1U)2| - \Wl\ 



MY \wi\ 



Dividing both sides of this expression by — ^ ' yields 



(kl' + KI 



2\ / 1 
1 + - 

2x 



2t;Jui ( 1 



2 2 

(1^21 + \vi 



1+7' 
4^2 



1 + 7' 

(hi' + Ki 



27^ 



(|W2|^ - |wi|^) , 



1 



47^ I |2 



v2 1 _,_ ^2 

1 + 7^ V 



7^ 
272 



2t;^i;i I 1 



1 + 7" ■ ■ i + 7 
Now the constant in this equation can be written as 



1 + 



l_ _ 2 \ _ 1 + 7^ 
72 1 + 7^^7^1 + 7')^ 
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SO that 



{\V2\' + \V,\ 

472 



1 + 7' 



7^1 + 7^) 

27' 



2vjvi ( 1 



1 + 7' 

Dividing by yields 



1-7^ 
1+7'' 
2 



(hl' + kii' 
47" I |2 , 27^ , , I 

and if 9"^ = then 1 - 6'^ = and so HA.2|) reduces to 

= (l-02)|^2K|V(|ii;2|'-hi| 
Some further rewriting shows that 



(1-1 



\V2+Vi\ + 



2 , (1 + ^') 



\V2 - Vi\' 



{l-e^)(2\w,\' + \w2f~\wif), 



or 



so that finally 





(1 + 




+ Wlp - 


^(l_02)l-2 ^ 


2 (2u;^ ^ 






^2 + Vi 


^ 1 




2 




2 



(A.2) 



since j^fr = This concludes the proof. 

Appendix B. 

Proof of Theorem H2.1() . This Theorem can be generalized slightly to the case 
where M and m are matrices. So if the collision matrix can be written as F = ( y ^ ) 
then the solution to equations H2.4a|l - H2.4c|) equal 

P ^X{m, M) + (Af + Mm-^M)~^ UY{ra, M)^ , 

G = X{M,m) + {m + mlVr^my^ UY{M,m)i, 
V ^ m-^M{I- P), 
Q = M-'^m{I-S), 
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with 

X{m, M) ={M + m)-^M, 
F(m, M) =M - Mm'^M 

+ X{m, M f {M + Mm-^M) X{m, M). 

This reduces to (|2.5|l if M and m become diagonal. Notice that solution H2.5|l can 
also be written as 

P = cos {ei2f {I - -f^U) , Q = sin {9/2f (I + U) , 



V = cos (d/2y {I + U),G^ sin {0/2y ~ —Uj , 
U'^U = I. 

To prove that (|2.5|l and (|B.1|) above are solutions to equations ()2.4a|) - (|2.4c|) 
notice that (|2.4a|) is equivalent to V — m^^M{I — P). Substituting that into 
P'^MP + V^mV = M yields 

P^MP +{I- P)'^Mm~^M{I - P) =M, 



or 



P^{M + Mm-^M)P - P'^Mm-^M 



so that finally 



(P - Xy^iM + Mm-^M){P - X) = Y, (B.2) 

with 

X{m, M) ={M + Mm-^My^Mm-^M = (to + M^^M, 
Y{m, M) - MnT^M 

+ X(m, Mf{M + Mm-^My^X{m, M). 
The solution to equation (|B.2|) equals 

P = X{m,M) + {M + Mm-^M)-iuY{m,M)^, (B.3) 
and V follows from V = m^^M{I — P) so that 

V = m-Hl (^I - X{m, M) - [M + Mm-^M)-^UY{m, M)^^ 
= m-^M {{m + M)-^m - {M + Mmr^M)-^UY{m, M)^^ , 
Specifically if M — MI,m = ml (i.e. the matrices m, M equal the constant 



(B.4) 



masses m, M times the unit matrix /) then (jB.2|) simplifies to 



27 



Then 

27 y V 27 
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Since 

.2/ 



1 ^ 

2 ^4 

cos{e) 



472 ' 4 

I V 7 



,1 + 7V (1 + 7^) 
Equation (|B.3ll then simplifies to 



for some unitary matrix U and so equation (|B.4p for the matrix V reduces to 

27 

Using the same approach the matrix G can be calculated as 

G = X{M,m) + (m + mA'r^myiUTY{M,m)i, (B.5) 

for another arbitrary unitary matrix U-\-. Again, if M = MI,m = ml (i.e. the 
matrices m, M equal the constant masses m, M times the unit matrix /), then 

7 sin(g) f 1_ 
2 7^ 



G=I^(/-4f), (B.6) 



„ ^sm(9) , ^ , , 

Q=^^(I + V). (B.7) 

Finally from the fourth equation 

= ^(1 + 7^)(/-C/-C/t), 

from which follows UUt = I so that U = Ut {U and C/t are unitary) and the proof 
is complete. 

Proof of Theorem 1)4. Using ||^J) and denoting A+w = + A^w = W2 —Wi 

it is clear that 

A-vJ \v2-vj \P-I QJ \wi 
P + I Q\ fvi 
-Q Q Ui 



so inverting yields 



vi \ _ fP + I Q\ ^ fA+v 



wi I \ —Q Q J \A V 



(B., 



which can be entered into the total kinetic energy expression Tik = vf+wf = V2+W2 
once the inverse of this matrix has been determined. Since P + Q = / it is clear 
that P + Q + / = 2/ so some manipulation shows that 



P + I Q\ ^ fl 



Q QJ 2\I (2Q-1 



(B.9) 
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Hence 



I I 

-I {2Q-^~I) 



I 

72 

/ 

7' 



P + I Q 
-Q Q 

I -I 
I (2Q-1-/) 



1 + 7' / 



where = {2Q-'^-l). 
As a result 



M 



7^ 



I 7^^!^ 



/ -A (A+v 
I Vt) \A-v 



A+v 
A^v 



(B.IO) 



/ (1 + 72) / -/ + j^Vr 

= (1 + 7^) A+v^A+v - 2A+v^ {I - -f^Vr) A'v 

+ A'v^ (/ + j^VjfVr) A'v 

To further simphfy the appearance of this expression the following Lemma is re- 
quired. 

Lemma B.l. Let Z be the collision scattering matrix Z = I ~ 2{I + U)~^ = 
I — 'ysm{6)Q~^ with U the unitary matrix and Q defined in equation IjB.Tp . Then 
it is true that 

1 + 7^ 



T 



Q 



{i + u)-\ 
2 



jsin{9) 



I, 



and 



Z'^ + z-^ = 0, 



= Z, 



1 - -f^Vr 

1 + 7' 
1 + 7'^r^^r ^ 



1+7-^ 7^ 
Proof. Using result (|2.5|l it is clear that 

2 



'-±Xz-z. 



(B.lla) 
(B.llb) 

(B.12a) 
(B.12b) 

(B.12c) 



U 



7 sm{9) 

for the unitary matrix U. Hence, 
I=U'U ■ ^ 



Q-I 



Q-I 



1 + 7^ 



Q 



7sin((?) 
4Q^Q 2 



7sin(6') 



Q-I 



')'^siv?{9) 7sin(0) 



(g^ + g)+/, 
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SO that 



2Q^Q (Q^ + Q) 



7sin(6') 

Multiplying left with matrix Q""^ and right with matrix Q^^ it follows that 

7sm(6') ^ ' 

which demonstrates equation ljB.llb|) . 

Substituting definition Q^^ = {I - ^)/ (7sin(0)) into (|B.llb|l yields 

2 _^_f{I-ZY , {I-Z) 



7sin(6') \ 7sin(6') 7sin(0) 



7sin(^^) 7sin(0) 

from which follow Z"^ + Z = proving Hij.l2a(l and equation (|B.12b|) follows from 

= (/ + 7^) (/-2(/ + [/)-!) 

= (/ + 7^) (/ - 7 sin(0)g-i) = (1 + 7^) Z. 

Finally to demonstrate IIB.12c|l use the definition Vr = 2Q^^ — I again to show that 

f I - Z \ 2 Z 

Vt^2{ -—- ] -1= . - / - 2- 



so that 



7sin(0) / 7sin(6') 7sin(6') 

/ + i^V^Vv = 1 + ^ (/ - (1 + 7^) Z^) (/-(! + 7') Z) 

which concludes the Lemma. □ 
Applying (|B.12a|l - HB.12c|l to (jKTUIl yields 

^.AW.-2A+.-(^^^A-. 
Mt (1 + 7^) 

+ ^ ^ (1 + 7^) ^ " 

=A+t;^A+t; - 2A+w^ZA"z) + ^tA^u^A^i; 

7 

A-z;^Z^^A-t;, 

which is equivalent to equation H4.1|l and this concludes the proof. 
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The total kinetic energy Hk can also be expressed in terms of the velocities of 
the heatbath particle and the same random matrix Z as is shown in the following 
Theorem. 

Theorem B.2. The combined energy of the main and colliding heatbath particle 
Tik can be expressed as 

^ ^A+w'^A+w - 2A+w'^ZA-w 

Mt (B.13) 

+ -f'^A^vFA-w + (l + 7^) A^vFz'^ZA-w, 

where 

A^W = W2 + Wl, 
A^W = W2 — Wi. 

As previously the anti- symmetric matrix Z is defined as Z = I — 2{I + U)^^ — 
I — J sm{9)Q^^ with U a unitary matrix and Q as defined in H2.5|l . 

Proof. Using (|4.1|l again it is clear that 

A+i(;\ _ /'w2 + wi 

A^W J \W2 — Wi 

V G + I\ fvi 

V G-IJ \wi 

V G + I\ (vi 

V -V ) \ wi 

so that 

vi\ _ (V G + (A+v 
wi) ~ \y ~V ) \A-v 

As before use y + G = / so that y + G + / = 2/ to show that the inverse can be 
calculated as follows 

V G + I 

V -V 

with X = {2V-^ - I) so that 

Y' fl 0\fV G + I 
G^ + I -V^J [o 77 \V ~V 

-/)(o 7°)(/ 

(1 + 7')/ X-72/ 
X^-^^I X^X + j^I 









)4 







(B.14) 
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Hence 



= ^ ™ ) (^>T „ ;^T;^ + j (^^-^ (B.15) 

= (1 + 7^) A+iy^A+w - 2A+w^ {-f^I - A^w; 

To further simplify the appearance of this expression the foUowing Lemma is re- 
quired. 

Lemma B.3. Using the same definition as in Lemma (jB.l|l let V^^ = "f'^Q^^ = 
^^^Tg. (/ + [/) = 7(/— Z) / (sin(0)) with U the unitary matrix and using equation 
ifi?.6)) it follows that 

^ {1 + j^) {I + Uy\ (B.16a) 

y-^ + y-i = -^/, (B.i6b) 

sin(d) 

and 

= Z, (B.16c) 



1 + 7^ 
X^X + 7^7 

1 + 7" 



= (1 + 7^) Z'^Z + 72/. (B.16d) 



Proof. Using V ^ = 7^Q ^, equations ljB.16a|) and (|B.16b|l foUow immediately 
from HB.lla|l and (|B.llb|l . Now X = 2V-^ - I = 2j^Q-^ - I, hence 

(X- 7'/) =2 (1 + 7^) (/ + C/)-i-(l + 7')/ 
= (1+7^) (2(7 + [/)-!-/) 
= -(1 + 7^) Z, 

so that X = 7^7 - (l + 7^) ^ which proves IjB.lGcp . Then finally 
X^X + 7^7 

= (7^7 - (1 + 7^) Z)^ (7^7 - (1 + 7^) Z) + 7^7 
= (l+72)'z^Z + 7'(l + 7')/, 
since Z'^ + Z = Q. This concludes the proof to the Lemma. □ 

Applying expressions (|B.16ap - HB.16d|) to ||5T5|) yields 

^ ^A+w'^A+w - 2A+w'^ZA-w 
Mt 

+ 72A"w^A"-u; + (1 + 7^) A^nF Z'^ ZA^w, 



which is equivalent to equation (|B.13|) and this concludes the proof of the Theorem. 

□ 
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Appendix C. 

Proof of Theorem 1)4.3(1 . This argument is a variation on the proofs presented by 
E. Nelson [7| and E. Carlen To prove Theorem (|2.12|l express the probabihty of 
the particle position as p = p{x,t),x € R" for an appropriate (smoothly) differen- 
tiable function R = R{x, t), x G M" such that p — e = e ^ . Then introduce the 
sufhciently smooth functions A = A{x, t), S ^ S{x, t), a; e M" and constants 6, ^ to 
express the backward and forward drifts fe+ = t), = b^{x, t),x G R", t > 
as follows 

2j5\'R^a^— = b+-b-, (C.la) 
P 

2^{WS - A) =b+ +b~, (C.lb) 



or equivalently 



6+ = ^ {VS -A) + -f6VR, 
6" = ^ (VS" -A)- -fSVR. 



(C.2) 



Equation l|C.la|l is a consequence of equation H1.18c(l and the fact that p — e 
while equation l|C.rb|) supplies a definition for the functions A = A{x,t),S — 
5(a;,t),a; e R". 

Using this definition equation 1(4. 4|l becomes equivalent to 



(C.3) 



2 



Mt 

2 



e I V5 - Ar - 2^(57 (VS" - Ay ZVR 
+6^VR^r^VR 



1 

2r ' 2x72 " ■ ' iir' 



where = /+ (1 + 7^)^^ [ZZ"^] and where Z is the collision scattering matrix such 
that E [Z] = Z. The time derivative of this functional depends on the continuity 
equation H1.19c|l which demands that 

Pt = -V. (^-^^p) = -eV. ((V5 -A)p), 



9/1 2~,SR 26E. 

where rj — a /"/ and p = e "'^ = e 1 . 



From this follows 



Rt — ^ {SxjXj — Ajxj) — i {Sxj — Aj) Rx- , 
Rtxk — {SxfXj ~ ^jajj) ~ i {Sxj — Aj) Rxj 



(C.4) 



using Einstein's convention for the summation of indices. This abbreviates Sx^xj 
^xS = (£| + ... + £|) and Ajx, := Ej = ^1^1 + •■•^na;„. 



ENERGY AND MOMENTUM CONSERVATION FOR STOCHASTIC PROCESSES 61 



Taking the time derivative of the kinetic part of the Hamiltonian H4.6(l using 
^ and ignoring the trace term results in 



d (\ 



dt \2 
1 



E 



' [s.,^t - i,) ^,fei?., + {Sx, - A,) ^,kR 



Pt 

£,Sj I p , . _ 
p( e{S.,-A,)[s^^^,~A,)+ \ 

so with rearranging this reduces to 



dx 



dx 



dx. 



dt\2 



M' 



■'IT 



e {Sx, - A,) (5,^. - A,) ~ 2i6j {S,^ - A,) Z.kRx 



-S^Rxj ^jk^Xk 



^Sj J p (^(^Sxjt - Aj^ ZjkRxu + {Sx, ~ Aj) ZjkRxk^ dx 

f e {Sx, - A,) {^S^^t - A,) 

3^R-:^jk - {Sx, - Aj) Z,k) Rt. 



\ 



\^^Rx,^%Rxk 



dx. 



dx (a) 
(b) 



The first term becomes 



(a) 



I +S'RxT^,kR^, 

'^(^2^ {^x:pXp - Apxj,) - i {Sxp - Ap) Rxp^^ 

(e{Sx,-A,){S,^_-A,)\ 
-2^S^ {S^^ - A,) Z.kRx, 
\ \ +5^RxT^,kRx. J I 

e {S,,-A,){S,^-A,) 
-2^5-1 {S,^ - A^) Z.kRx, I I dx 



dx 



6 

= - I P 



. ■ ~ A 1 



~^^Rxj^%Rxk 



.5 r ( ({e{Sx,-A,){S,^^A,y 
-— P\{K- A) -2e<57 {S^, - A,) ZjkRx 

^ \ V +S^RxT^,,Rx, 

so with one partial integral this reduces to 

/ f (Sx, - Aj) {S^^ - Aj) 

(a) =Up\ {Sx, - Ap) -2^Sj\s.^ - A,)%kRx, 
V V +S'RxT^,,Rx, 



dx. 



dx 



I / - ^p) -2^'57 {Sx, - A,) ZjkRx, dx 

^ ^ ^5^RxT^,,Rx, II 
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-- P\{S.,- Ap) R., -2C<57 (5., - A,) Z,kR., I I dx, 



or 



+ ^ / -pR^ 

2 J J? 



da; 



\ -J jn. / jjp/ 

{Sx, - Ap) -2^(57 (^x, - ji) Z,kR., 
\ \ +6'^R^jr^-^Rx^ 

-- P\iS-.- A) -2?^7 (^x, - A,) Z,kR 



SO that finally 



dx 



dx, 



( e{S,,-A,){S,^-A,) 



Now the third term becomes 



dx. 



e (.Sx, - A,) (S^^t - Aj) + \5^R^V 

(c) = / p I V <5'i?x,r|, _ ^^ ^-2£ (5^^^^ _ A 

-^5^[S,^-A^)Z, 




'jTs/ — Aj^ R, 

( e (Sx, - Aj) (5,,t - i,) + i^^i?,, f 



— Mr 



-C<57 {Sxj - Aj) Zjk 



25 ("^^p-p 

C ('S'^p ~ ^p) Rx 




— A 



= p 



^ ■e<57 (5x, - A^) ZjkJ V (^xp - Ap) R,^ J 

e{Sx,-A,) (^S^^,-Aj) + '^S^R,T^ " 



+ if (<52i?xT|, - ^<57 (^x, - A,) Z,,) (5.^.^ - A 

\ / Xfc 

H (s'Rx, r|, - ^,57 {S,^ - A,) Z,A {S,, - Ap) R,^ 

HRx, [S'Rx T^.k - (^., - A,) Z,k) (^xpxp - 
-'-^ (<52i?, T|, - ^(57 (^., - A,) {S,^ - Ap) R,^ 



^pxp, 



da;. 



This can be written as 

/ e {Sx, - A,) (5,^., - i,) + i(52^,^ f \ 

(c) = y p +e (-J'i^x, r|, - e<57 (^x, - A,) Z,k) (^xp - Ap) R^^ 

^ ^{^'R^,^h-^h{Sx,-Aj)-Zjk)^^{Sx,x,-Ap,^) \ 



(c) = / P 



.dx 
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and with one partial integral again this becomes 

H (s^R., r\ - ^Sj (S,^ - A,) Z,k) (S., - Ap) R,^ 

, -^^R'^. [s^R-,^lk - {s., - A,) z,u) (5., - Ap) R,^^ 

( -^R,Jd^R.,T%-^d^{S^.-A^)Zik) {S.,-Aj) \ 
P -?^x.x, {6^R., r|fe - (5x, - A^) Zjk) (5., - Ap) 
^-^^R,^R,^ (^'i?x,r|, - (5., - Aj) Zjk) (S., - A. 

f /-il('5'^x,r|,-^^7(5x,-^,-)^ifc) 

So then combining the a), 6) and c) terms results in 
(a) + (b) + (c) 
/ 



+ 



+ 



dx 



■Ap)J 



dx. 



e (5., - A,) - 4) + ^S^R,^ t^^.R,, \ 
'^Rrr.T'. - £<57 (5,, - A,) Z,^.) (S^^ - A„) 



+ 



-'V-./ ^/\-J- J/ ^ —jjri.-ft. 

('^'^-.r|fe - ^<57 (5., - A,) Z,k)^^ (s., - Ap) 
y+^r^.. {s'RxT^jk - {Sx, - A,) Z,k) (5., - Ap) R,^^ 
I -^R^J5^R.,T^.^-^5^{S.,-Ai)Zik) {S^^-A^) \ 
-Ci^x.x, {^^R.,T% - {S,^ - A^) Zjk) {S,, - Ap) 
-^-fR^^.R., (-^'^x.ri, - ^^7 - A^) Zjk) (5x, - A,)^ 

f (-^(^'R^.^]k-^^iiSx,-A,)z,,) 

IP dx 

J y-S,R., (^'^x - ^(57 (5., - A,) {S,^ - Ap) j 

-^h j P ((-Sx,* - Aj) ZjkRx^ + {S^j - Aj) tjkRx,) dx 

, ,co / /e{S.,^A,){S,^^A,)' 
7i CP (^x, - Ap) -2^(57 {S,^ - A,) ZjkR., 
^J-oo \ \ +5^R^T'.^R^. 



dx 



dx 



+ 

which reduces to 

(a) + (b) + (c) 

p(e {s.,-A 



J P {e {S., - A,) (5,^.* - Aj) + ^S'R^T^^.R,,) dx 
-^jp{S.,- Ap) {R,, {S^R,,T^jk - (5x, - Aj) Zjk))^^ dx 
-%j P{S.,- Ap) (<5'i?x,q, - ^<57 {S,^ - A,) Z,,)^^^^ dx 
-^h j P {{s^it - ij) ZjkRxk + [Sx^ - Aj) tjkRxk^ dx 
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p {S,^ - ^p) ( -2^57 (5,^. - A,) Z,kRx, 
+5^RxT^,kRx, 



dx. 



Gathering terms finally shows that 



dt 2 

< I P {Sx^ - Apj 



) 

iSt + 3" {S^. - Aj^ {S,j-. - Aj 



Srj 



dx 



J P (.Sx, -Ap)Ap + -J p [S^R,T^^,R,,j dx 
+ 1 / P (^.. - ^p) ((^^7 {Sx, A,) Z,,),,, J dx 
-i5^ j P [(S^x.t - i,) Z.kRx, + {Sx, - A,) t.kRx,) 



dx. 



(C.5) 



Now the static potential term time derivative j^-^E [^p] contains the differ- 



ential -n—E 



{\IS~Af A 



The assumption on the potential turns it into a full 

differential and the A term disappears if ^ = -^f^ ■ Once the term with the Trace of 
the matrix is reintroduced the Theorem is proved. 



Therefore a sufficient condition for the Hamiltonian in proposition to be- 
come time-independent the integrand in 1)4. 6|l is to equal zero. The following propo- 
sition shows that this is equivalent to demanding that the probability for the main 
particle position p{x, t) is derived from a wave function satisfying Schrodinger's 
equation. 

Proposition C.l. Assume that Z = so that E [ZZ'^] = 0, E[Z] = 0, = / 
and let the potential $p satisfy 1)4. 5|l . Define the wave function 

4> — ip[x, t) = e ^ , 
with p{x,t) = \%l!{x,t) p. Then equation 

eiSt - At) - ? |Vi?|2 + ? - A\^ - ^A,i? + = 0, (C.6) 



2 ' '2 (,AL 



T 

makes the energy terms invariant. Here Ax — \ '§^^ ) ■ Now this equation 

equivalent to 

1 2 
«xV'i = -XTi- (xV - iA) ^p + (j){x,t)^p. 

with X — Mxri and S = ^ = \/Mt so that ip — il^{x, t) — e x 
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Proof. The proof involves a straightforward verification of equation 1)4.11(1 by brute 
force. Taking derivatives 

I 

Tj 1 

Tj^ 1 2 ^ 

i^r;^,Vx^ = ip {-iS^A.Rx^ + A.^^Sx^) , 
which combined becomes 

= (^SRt - 5^A,Rx^ + SiRx^Sx, - "^Ajx, + y5,,^ 



Now 



= - ^ {v^i^x, - rjiAjS^ip) + ]^iAj^ {rjipx^ - iAj^-ip) 

= -\ W"^x,x, - ^iAjiijx, - iii^Ajx^ip) + ^iAjS,'q^ljx^ + ^ \Af fi}} 

= - ^v'^i^x.x, + ^iS,Ajx^i; + iAjS^-qipx^ + ^ |v4|^^^V, 
so indeed 

iV^t = (»7V ~ iAS^f tp + j^-fp- 
Finally, to adjust for the masses, let x — M^t] and let ^ — l/M^, then 



and 



R(x,f) + »S(x,t) 

= yj{x, t) = e X 



1 2 



This proves (|C.4|I . Also quite clearly 

b+{x,t) = ^i\/S-A + ^Rx), 

b-{x,t) = ^(VS-A-jRx), 

which proves equation H2.18|l . □ 

An immediate conclusion from quantum mechanics is that the mean acceleration 
of the motion is guided by the potential </>. 
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Proposition C.2. Some straightforward derivatives show that 

^^EHt)]--^Jpm)d.. 

Proof. In addition, 

j^E[x{t)] =-i j xV. ((V5 - A)p) dx = i j pV{S - A)da 



while 

72 



^E[x{t)] j pV{S ~ A). {A.^S - V.^) dx + i j p{VS - A)dx 
^ J (^^St + ^\{VS~A)\'^ dx 

p (V0) dx. 



1 



□ 



Now the example (|4.1|l can be slightly extended in the form of the following 
Proposition. 

Proposition C.3. Consider the two-dimensional case n=2 and assume that 

so thatE[ZZ'^] = 1/2/ andV' = I+{l+-f^)E [ZZ'^] = (l + [l + I = all. 
Then define the wave function 

^fl(x,t) + iSS(x,f) 

tp = ip^x, t) = e 1 , 
with p(x,i) = \jjj{x,t) p. Then equation 

makes the energy terms invariant. Here A.^ = {^-§^i ^3"^- -^ow this equation is 
equivalent to 

ix^H = -^{x^-tAf,p + <Pix,t)i, (C.8) 
with X = Mtii/<Ju, S, = 1/Mt and S =1/ {gj^Mt). 

Proof. Let the wave function 'ip = ^'(2;, <5, x) satisfy equation (|C.6p . Then it is 
clear that the wave function -0' — ip{x,t, — — ) satisfies equation l|('.7(l above. 

□ 
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Appendix D. 

Proof of Proposition (14.5(1 . Move all the main particle terms in equation 1(2. 2|l 
to the left and all heatbath particle terms to the righthand side (reference also 
equation 1(2.5(1 ) to obtain 



(D.l) 



/ -P\ {V2\ ^{ Q\ (W2 

-R) \yi) \~I S) \wi 
and using remark ((2.2(1 this can be reduced to 

vi) [O -R 

I -PR-^ 


^j'Q-^ 

iHQ-'-i) -i'Q-' + {i + i')i\ fuj2 

Using the collision scattering matrix {I — Z) — 2jj^^^Q^^ the expression above 
reduces to 




V2\ _ 7 /cos(6') 1 \ (W2 
vij ~sm{e) V 1 cos(6l)y \wi 

7 f-Z Z\ fw2 



s\n{9) \-Z Z) \wi 

V2\ _ 7 /COS(^) 1 \ (W2\ ^ fW 



vij sm{9) V 1 cos(6')y \wi J \W 



(D.2) 



(D.3) 



where — mW — '^t^ mZ{w2 — wi) — 'mZA~w is the additional momentum 
transfer from the heatbath to the main particle. Notice the rather interesting fact 
that this allows ((D.3(l to be rewritten as follows 

'V2\ ^ {C0s{e) 1 \{W2+ W\ . . 

^vj sm{9)\ 1 cos{9)J \wi+WJ ' ^ ' 

since the vector ( ]^ ) is an eigenvector of the matrix ^j^^ cos{e) ^J^^^^ ^ This is in 
fact equivalent to equation ((4.1t)() . 
Now returning to ((D.2(l 



V2\ _ 7 fcos{9)I -Z I + Z \ fw2 



vij sin{9) \ I-Z coa{9)I + Z J \wi 



so that 



V2 



W2\ _ sm{9) fcoB{9)I-Z I + Z 
wi/ 7 1^ I-Z cos{9)I + Zj \vi 

1 /- cos(6l)/ - Z I+Z \ fv2 



7sin(6l) V -cos{9)I + ZJ \vi 
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Then 



E 




Tl 



2a2 



E 



cos{d)I -Z I + Z 

i-z -cos{e)i + z 

- cos{0)I -Z^ I- Z^' 

/ + - cos(0)/ + Z'T 

(1 + cos2(6i)) / -2cos{9)I 



r7sin (6) 



-2cos{e)I (1 + c 

E [zz^] n 

E [zz^] 



where 



n 



(1 - cos{0)I) Z + E [ZZ^] . 

Recombining the cos{9)^ and sin(6')2 terms in the first matrix (refer to Theorem 
(12. SI) ) this reduces to 



E 





















;*) 


V T2 













ith 



2a2) 



(1 - 2a-] 
1 



2a' 



E [zz^] n 
n'^ E [zz^] 



r7 sin^ {9) 

Clearly here T z is positive definite and symmetric. This concludes the proof. 
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